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Resumo 



Investigamos o surgimento do comportamento coletivo global em redes de osciladores idênticos 
difusivamente acoplados, que no modelo estabelecido é uma variedade invariante pelas equações do 
movimento. A inter ação é modelada através da teoria de grafos e de sistemas dinâmicos. Utilizamos 
a teoria de contrações uniformes em equações diferenciais lineares não-autônomas para estabelecer 
os critérios sobre o parâmetro global de acoplamento que por sua vez define o estado síncrono 
e sua respectiva estabilidade sob perturbações lineares e não-lineares. O parâmetro crítico global 
de interação é dado em função somente da dinâmica individual dos osciladores, de propriedades 
espectrais da função de acoplamento e do segundo autovalor do laplaciano da rede. 
Palavras-chave: Sincronização, Redes, Estabilidade. 
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Abstract 



We investigate emergence of the global collective behavior in networks of diffusively coupled 
identical oscillators, which in the established model is an invariant manifold of the motion equations. 
The interaction is modeled with the graph theory and dynamical systems theory. Wc use the uniform 
contractions theory in non-autonomous linear differential equations to address the criterion under 
the global coupling parameter, which it turns defines the synchronized motion and it stability under 
small linear and non-linear perturbations. The criticai global interaction parameter is given only by 
the isolated dynamics, by the spectral properties of the coupling function and the second eigenvalue 
network laplacian. 

Keywords: Synchronization, Complex networks, Stability. 
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Capítulo 1 

Introdução 



o termo sincronização, que do grego significa "ocorrência ao mesmo tempo", está relacionado a 
uma gama de fenómenos presentes em muitos ramos das ciências naturais, engenliaria e vida social 
[APK01],[Str03],[PBK07]. Este fenómeno está enraizado, por exemplo, na vida humana, ocorrendo 
desde processos metabólicos em nossas células às mais altas tarefas cognitivas [ADGK^OS]. Um 
dos primeiros indivíduos a estudar esse tipo de fenómeno foi o físico alemão Christiaan Huygens 
(1629 - 1694). Ele descobriu que relógios de pêndulo, quando pendurados em um mesmo suporte, 
tendem à sincronização (movimento dos pêndulos). De forma geral, esse estado síncrono surge da 
competição e colaboração entre os elementos de uma rede. Exemplos vão desde a sincronização em 
relógios de pêndulos à doenças neurais como Mal de Parkinson [TRW^98] e Epilepsia, onde está 
última acontece quando um grupo específico de neurónios mantém um estado síncrono [MJ03]. 

Modelamos então essa interação entre elementos através de uma rede, em termos matemáticos, 
através da teoria dos grafos e da teoria de sistemas dinâmicos. Passamos a distinguir a estrutura 
da rede, a natureza de interação entre os elementos e o comportamento dinâmico individual dos 
osciladores (elementos da rede) . Uma rede é dita complexa quando a mesma não possui uma estru- 
tura regular de conectividade, e no mundo real a análise de sistemas interagentes é feita a partir 
da modelagem em redes com estrutura complexa. Existem vários tipos de redes que se enquadram 
nesta característica, como por exemplo as redes aleatórias, redes com estrutura de pequeno mundo, 
redes scale-free, entre outras [CLotMSU06]. 

Estaremos interessados em sincronização que ocorre de forma global. E global no que se refere 
a totalidade dos elementos da rede e ao tempo, ou seja, uma vez que o estado síncrono ocorre, 
o mesmo permanece por todo o tempo futuro. Um dos principais ingredientes para o surgimento 
da sincronização global na rede é que os sistemas interagentes sejam idênticos, ou seja, que todos 
eles sejam descritos exatamente pela mesma equação de movimento [RP03]. Estaremos interessados 
apenas em redes de sistemas idênticos. 

Recentemente, Pecora e Carroll [PC98] estudando modelos de osciladores em redes comple- 
xas com interação semelhante â difusão, demonstraram que tais redes de sistemas dinâmicos não- 
lineares, porém idênticos, podem sincronizar globalmente, mesmo manifestando uma dinâmica in- 
dividual complicada. Utilizamos um modelo de rede semelhante ao utilizado por Pecora e Carroll 



Considere que um grafo G modela uma rede de n vértices. Sobre cada vértice da grafo intro- 
duzimos uma cópia da equação diferencial x = f{x), f : M™ — >• M"*, onde requeremos que / G C"', 
r > 2, de forma que, essencialmente, / representa uma dinâmica não-linear, possivelmente caótica, 
e que as condições de existência e unicidade sejam satisfeitas. Para o acoplamento entre os osci- 
ladores, consideramos o caso de interação difusiva, isto é, o acoplamento depende da diferença de 
uma função dos estados dos vértices que estão conectados. Então, dado um vértice i em G, a sua 
dinâmica passa a ser descrita na forma 



em [PC98]. 



n 
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onde os Aij^s são as entradas da matriz de adjacência, que codifica a informação topológica da 
rede dizendo quais elementos possuem interação direta, a é o parâmetro global de acoplamento e 
H : M"* — )• é uma função de acoplamento, que sem perder a generalidade podemos considerar 
que H, é uma matriz. Trabalharemos apenas com o caso em que H é uma matriz positiva-definida. 

Note que um importante detalhe deste modelo é que quando todo i,j, então 

o termo de acoplamento desaparece identicamente de forma que Xi{t) = Xj{t) ficam iguas por 
todo o tempo, caracterizando assim um subespaço invariante. Passamos então a chamar o estado 
síncrono de variedade de sincronização. A variedade de sincronização é uma variedade invariante 
pelas equações do movimento e a análise da sincronização global mostra que tal movimento síncrono 
e sua consequente estabilidade é definida pela intensidade do acoplamento a. Objetivamos então 
caracterizar a estabilidade desse subespaço, que possui uma geometria trivial. 

Nos últimos 20 anos houve uma evolução significativa no campo da sincronização em redes. A 
análise da sincronização beneficia-se da análise estrutural de tais redes [BS09], isto é, a estrutura 
de interação entre os elementos influência significativamente na sincronização [PC98],[Wu07]. Por 
exemplo, sabe-se que redes com estrutura de pequeno mundo, ou seja, redes bem conectadas, pos- 
suem uma maior propensão à sincronização do que as redes regulares que as geram, e consequente 
menos propensas à sincronização do que as redes complexas aleatórias advindas da construção da 
rede pequeno mundo [WS98]. Por outro lado, a heterogeneidade na rede, isto é, redes em que a 
maioria dos elementos possuem poucos vizinhos e apenas alguns elementos com muitos vizinhos, 
dificulta a sincronização global da rede [MZK05], no entanto facilita a sincronização entre elementos 
que possuem uma grande quantidade de vizinhos [PcrlO]. 

O modelo (1.1) tem sido muito usado como ponto de partida para se entender a emergência da 
coletividade. Nijmeijer e Pogromsky em [NijOl] estudaram o modelo do ponto de vista de controle 
e eles provam que as soluções existem globalmente. Hasler e colaboradores em [BBH04] estudam o 
problema construindo funções de Lyapunov para provar a estabilidade da variedade de sincronização. 
Macau e colaboradores em [GMY12] também estudam esse problema utilizando a existência funções 
de Lyapunov para garantir a estabilidade do comportamento síncrono. Josic em [JosOO] estudou esse 
problema do ponto de vista das variedades invariantes. Pereira em [Pcrl2] estudou este problema 
do ponto de vista de dicotomias e nossa abordagem é baseada na mesma perspectiva. 

1.1 Principais Resultados 

Os estudos aqui desenvolvidos provém uma compreensão rigorosa do comportamento coletivo 
em redes e assim abre frentes para novos resultados na área. 

Uma abordagem muito utilizada para se estudar a estabilidade da variedade de sincronização em 
redes é baseada na teoria dos expoentes de Lyapunov [PC98], [HCLP09], [PcrlO]. No entanto, sabe- 
se que neste contexto, perturbações arbitrariamente pequenas podem causar mudanças abruptas na 
estabilidade do movimento síncrono [BP02]. Abordamos então o fenómeno da sincronização através 
de uma nova perspectiva, isto é, através da teoria de contrações uniformes em equações diferenciais 
lineares não- autónomas, de forma que poderemos garantir que a variedade de sincronização seja re- 
sistente à pequenas perturbações. Destacamos as principais contribuições abordadas neste trabalho: 

1. Existência global das soluções de (1.1): Construímos uma função de Lyapunov para a di- 
nâmica isolada e extendemos para a dinâmica coletiva. Dessa forma, independentemente da rede 
considerada, as soluções de (1.1) são limitas e portanto existem globalmente. Veja o Teorema 16. 

2. Estabilidade da variedade de sincronização: A partir da teoria de contrações uniformes em 
equações diferenciais lineares não-autónomas poderemos determinar um parâmetro crítico de aco- 
plamento ctc tal que para a > ac em (1.1) a variedade de sincronização é estável. Precisamente, 
temos o Teorema 17. 

3. Persistência da sincronização: Introduzimos um modelo de perturbação que age na função 
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de acoplamento de (1.1) e utilizamos um resultado de persistência de [Cop78] para estabelecer a 
magnitude das perturbações que não destrói a estabilidade da variedade de sincronização. Veja o 
Teorema 18. 

Agora, gostaríamos de expor algumas considerações sobre os itens acima. No item 1 contruí- 
remos uma função de Lyapunov para a coletividade supondo que existe uma função de Lyapunov 
quadrática para a dinâmica isolada. O item 2 é um resultado fundamental e vários autores já o pro- 
varam [JosOO], [NijOl], [Pcrl2]. A novidade então, é que abordamos na prova a teoria de contrações 
uniformes em equações diferenciais lineares não-autônomas [Cop78]. Essa teoria não é usualmente 
empregada para estudar redes, porém foi recentemente usada para provar o item 2 em [Pcrl2]. 
O que fazemos aqui é expandir e abordar com mais detalhes a prova. Portanto, não há nenhuma 
técnica nova mas no contexto de sincronização o resultado é novo. 

O item 3 é o ponto novo e segue praticamente como um corolário do resultado apresentado no 
item 2. O caso onde / em (1.1) não é igual foi tratado por Pereira e colaboradores em [PEB"''13], 
porém o caso em que H não é igual é novo na literatura. No caso da pertuação na função de aco- 
plamento, a estrutura da rede tem um papel muito importante para a conservação da estabilidade. 

1.2 Organização do Trabalho 

Queremos deixar claro ao leitor que, do Capítulo 2 ao Capítulo 5 temos por objetivo fundamentar 
toda a teoria matemática que será utilizada para modelar o nosso problema, a saber, o problema de 
estudar a o surgimento da sincronização e sua estabilidade em redes. Portanto, o leitor deve saber 
que esse é um tema interdisciplinar e que por essa característica, faz-se necessário fundamentar 
e então interligar algumas diferentes áreas da matemática, para então dar o devido tratamento 
rigoroso sobre o mesmo. Isso faz com que, a princípio, pareça que estamos tratando de temas 
que estão a parte do título apresentado pela dissertação, porém os mesmos são de fundamental 
importância para o tratamento do problema proposto. 

No Capítulo 2, apresentamos os conceitos preliminares, e as notações matemáticas que utili- 
zaremos durante o desenvolvimento do trabalho. No Capítulo 3 abordaremos alguns tópicos da 
teoria qualitativa de equações diferenciais ordinárias, à qual é usada para modelar e entender o 
comportamento individual dos elementos da rede. Abordaremos em especial a construção da função 
de Lyapunov. Em seguida, no Capítulo 4 trabalharemos sobre a estabilidade e suas persistências, 
da solução trivial em equações diferenciais lineares não-autônomas. No Capítulo 5 abordaremos os 
fundamentos da teoria espectral dos grafos, que por sua vez é utilizada para modelar topologica- 
mente uma rede qualquer. Seguido pelo Capítulo 6 o qual introduziremos formalmente o modelo de 
acoplamento difusivo em redes, e enunciamos os principais resultados da presente dissertação. No 
Capítulo 7 trabalharemos vários exemplos para a teoria abordada, principalmente com o exemplo 
mais trivial de rede, a saber, o de dois osciladores difusivamente acoplados, do qual extrairemos 
valorosos resultados. E finalmente, no Capítulo 8 apresentaremos as provas dos principais resultados 
enunciados no Capítulo 6. Por fim, no Capítulo 9 discutiremos algumas conclusões e abordaremos 
ideias sobre uma possível continuação do respectivo trabalho. Nos Apêndices A e B, trataremos 
de resultados gerais e específicos sobre respectivamente Álgebra Linear e Equações Diferenciais 
Ordinárias, os quais são utilizados ao longo da dissertação. 



Capítulo 2 

Conceitos Preliminares 



o presente capítulo é utilizado, principalmente, para introduzir a notação que servirá como 
base para todo o texto. Estão disponíveis nos Apêndices A e B outras informações basilares para 
o entendimento da problemática exposta. Trabalharemos essencialmente com espaços que possuem 
norma e produto interno e matrizes quadradas. 

2.1 Espaços Normados e Espaços com Produto Interno 

Considere M o corpo dos número reais. Indicamos por M+ o conjunto [O, +00). Seja E um espaço 
vetorial sobre M. Uma norma em E éa função || • || :£'—)■ M+ satisfazendo as seguintes propriedades: 

(i) Para todo x £ E \ {0} tem-se que ||x|| > 0. 

(ii) 11x11 = O se e somente se x = O (vetor nulo). 

(iii) Dados a £ M e X £ E, tem-se que \\ax\\ = |a|||x||, onde \a\ é o valor absoluto de a. 

(iv) (Desigualdade Triangular) Para todo x,y E vale que ||x + y\\ < ||x|| + \\y\\. 

Dado X = (xi, • • • , Xn) G M", importantes normas que surgem são 

1. (p-norma) ||x||p = (X^^Li l^^íl'')^''^; Quando p = 2 temos a norma Euclidiana. 

2. (Norma do máximo) ||x||oo = max, |xj|; 

3. Norma da soma {p = 1). 

Dadas duas normas || * e || • \\y em E. Dizemos que tais normas sao equivalentes se existem 
números positivos ci e C2 tais que 

Cl||x|||^ ^ — C2||^||ii) 

para todo x £ E. 

Devido a forte garantia do teorema a seguir, por todo o presente texto, não fazemos uso de uma 
norma específica na maioria das situações em que trabalhamos com tal conceito. 

Teorema 1. Duas normas quaisquer em M" são equivalentes. 

Prova: Ver [LimlO], p. 19. 

De forma mais geral, em espaços de dimensão finita, todas as normas são equivalentes [O'S06]. 

Uma norma em um espaço vetorial E dá origem à noção de distância em E. Para todo x,y £ E, 
a distância de x a y é definida por 

d{x,y) = \\x - y\\. 
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Dessa forma, como em um espaço vetorial normado tem-se a noção de distância, ou métrica, o 
mesmo pode ser chamado de espaço métrico. Um espaço métrico E é dito ser completo se toda 
sequência de Cauchy em E é convergente em E [O'S06]. 

Definição 1. Seja E um espaço métrico completo e f : U C M — )■ uma função contínua. Dizemos 
que f é uniformemente limitada em t, t £ U, se existe alguma constante c tal que 

sup||/(í)|| =c. 

íel/ 

Seja E um espaço vetorial sobre o corpo M. O produto interno é a função real (•, •) : ExE ^ M, 
que associa cada par de vetores u,v € E ao número (n, v), de modo que, são válidas as propriedades 
de simetria, bilinearidade e positividade [LimlO]. 

2.2 Matrizes 

Indicamos por Mat(M, n) o conjunto das matrizes quadradas de dimensão n x n com entradas 
em M. 

Definição 2. Dizemos que uma matriz H G Mat(K,n) é simétrica se H = H\ onde f representa 
o transposto. 

A matriz identidade, indicada por € Mat(M, n), é uma matriz simétrica. Qualquer matriz 
A G Mat(M, n) pode ser decomposta na forma 

, A + A^ A-A^ 

A = \ , 

2 2' 

onde a parcela [A + A^) /2 é chamada de parte simétrica de A. 

Reservamos a notação "diag" para representar as matrizes diagonais em Mat(M, n). Ou seja, as 
matrizes cujos os valores fora da diagonal são nulos. Por exemplo, I3 = diag(l, 1, 1). 

Definição 3. Seja x G \ {0}. Dizemos que uma matriz A £ Mat(M, n), simétrica, é positiva- 
definida [positiva semi-definida] se {x, Ax) > O [{x, Ax) > 0] e negativa-definida [negativa semi- 
definida] se {x,Ax) < O [{x,Ax) < 0]. 

Note que (x, Ax) = x^Ax, de forma que em muitos contextos poderemos usar essa equivalência. 
2.2.1 Teoria da Representação 

Definição 4. Diremos que uma matriz Q G Mat(M, n) é ortogonal se = Q^^, onde Q^^ é a 
inversa de Q. 

Definição 5. Uma matriz real A é ortogonalmente diagonalizável se existe uma matriz ortogonal 
Q tal que Q^AQ = D é uma matriz diagonal. 

Teorema 2 (Teorema Espectral). Uma matriz A G Mat(M,n) é ortogonalmente diagonalizável se 
e somente se é simétrica. 

Prova: Ver [Lim09], página 167. 

Corolário 1. Seja A uma matriz real e simétrica. Então existe uma matriz ortogonal Q tal que 

1. D = Q^^AQ é uma matriz diagonal onde os elementos da diagonal são os autovalores de A e 

2. As colunas de Q são formadas pelos autovetores de A. 

Proposição 1. Seja A G Mat(M, n) uma matriz simétrica. Então para todo X, autovalor de A, tem- 
se que A é positiva- definida [positiva semi-definida] se e somente se A > O [A > 0] e negativa-definida 
[negativa semi-definida] se e somente se A < O [A < 0]. 
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Como A é simétrica então podemos formar uma base de M" com os autovetores de A os quais 
podem ser tomados ortonormalizados. Seja /3 = {vi, ■ ■ ■ , u„} tal base, então x = Y27=i '-í^*' ^ 
o que produz 



{x,Ax) = (y2ciVi,Ay2ciVij 

= ( ^ CiVi, ^ CiXiVi ) = ^ Cí^Aí 



\ 2 l 11 

Sendo A» > O para todo i tem-se que (x, Ax) > 0. Os demais resultados seguem alterando-se as 
hipóteses sobre o sinal de cada Aj. 

2.2.2 Norma Induzida de Matriz 

Se II • II é uma norma em M" e j4 G Mat(M, n), definimos a norma de A por 

ll^ll = sup ||^a;||. 

xÇM.'^,\\x\\<l 

Essa norma é conhecida como norma induzida [LT85]. Por toda a dissertação, fazemos uso dessa 
norma ao se trabalhar com matrizes. Uma notação mais precisa seria ||^||p, onde o índice p indica 
qual a norma sobre estamos induzindo a norma em Mat(M, n). Omitiremos tal índice, indicando-o 
apenas quando for necessário. 

Observação: a noção de função uniformemente limitada dada pela Definição 1 pode também ser 
extendida à matrizes considerando a norma induzida. 

Uma das propriedades mais importantes dessa norma é a propriedade sub-multiplicativa. Mais 
especificamente temos o 

Lema 1. Dados A, B ^ Mat(M, n) temos 

1. \\Ax\\ < ||^||||a;||, para cada x £ M"; 

2. II Ai?|| < ll^ll ||i?|| (propriedade sub-multiplicativa) . 
Prova: Ver [DL07], p. 351. 

Exemplo 1. Sejam A G Mat(M, n), A = [Aij], e x £ M". Considere a norma do máximo definida 
em M" por \\x\\oo = maxj |xí|, onde x = (xi, • • • , x„). Então 

ll^lloo ~ sup ||^x||oo 

II II oo — 1 

^ max ^ ^ A^/i j X j 
j 

< max^^ l^jjl l^j I < ||a;||oo max^"^ | Aj| = max^^ |^jj|. 
j j j 

Por outro lado, podemos mostrar também que ||A||oo ^ max, | Aj | • Certamente, para algum 
1 < k < n tem-se maxj Ylj l^ijl — Ylj l^fcjl- Então, existe xq = («i, • • • , a„) G C"' com ||xo||oo = 1 
escrito da forma 

í ^ se Akj + O 
I O se A^j = O 
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que cumpre 



Portanto, 



|^a;o||oo = inax 



EA 



max 

i 



1^ 



kj\ 



A 



kj 



\A\\ 



max 



'^\Akj\ >max'^\Aij\. 



ax ^ ^ I Aij I , 



que simplesmente é a soma máxima dos valores absolutos das linhas de A. 

2.2.3 Produto de Kronecker 

Sejam A = [Aij] G Mat(M, n) e B = [Bij] € Mat(M,m), então o produto de Kronecker de ^ e 
B, escrito como A B £ Mat(M, n x m), é a matriz formada pelas submatrizes AijB escrita da 
forma 

'AiiB A12B ■■■ AmB' 

A21B A22B ••• A2nB 

A0B = 

AnlB Am2B ■ ■ ■ AnnB 

Teorema 3. Sejam A G Mat(]R,n), B £ Mat(M,m), U G Mat(M,n) eV e Mat(M,m). Então 

{A ^B){U(g)V) = AU(S) BV. 

Prova: Ver [LT85], p. 408. 

Teorema 4. Sejam A £ Mat(M,n) e B e Mat{R,m), então 

{A0 B)^ = A^ bK 
Prova: Por definição tem-se {A (g) B)^ = (Aj = AjilB^ = A'^ B\ ■ 
Teorema 5. Sejam A G Mat(]R,n) e B £ Mat(]R,m). 

1. Se A e B são positiva- definidas então A(^ B é positiva- definida. 

2. Se A é positiva- definida e B é positiva semi-definida então A(^ B é positiva semi-definida. 

3. Se A é positiva semi-definida e B é positiva definida então A0 B é positiva semi-definida. 

Prova: Considere o item 1. Sejam x G e y G não-nulos. Então por hipótese tem-se que 
x^^Ax > O e y^By > 0. Assim x^Ax ^ y^By = (x (8> y)'^{A ^ B)[x y) > O (Teoremas 3 e 4). Os 
demais itens seguem de forma análoga. ■ 



Teorema 6. Sejam x G M", y G 



\p uma p-norma. Então 



y\\p — ll^llplly llp- 

Prova: ||x(g)y||^ = ELi E/=i kfcy/T = ELi kfc|^E/=i IvÚ^ = 

O resultado é estensível para matrizes condiderando a norma induzida. 
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2.3 Conjuntos 

A bola aberta de centro em xq € M" e raio ô > O é o conjunto Bs{xo) = {x £ M" : d{x, xq) < 6}. 
Analogamente, a bola fechada de centro xq e raio 5 > O é o conjunto Bs{xo) = {x G M" : d{x, xq) < 
5}. 

Definição 6. Dizemos que um conjunto E C é limitado quando o mesmo está inteiramente 
contido em alguma bola fechada. 

Dado um ponto a £ E, dizemos que a é ponto interior a E quando, para algum 6 > O, tem-se 
que Bs{a) C E. 

Definição 7. Um conjunto E é aberto quando todo ponto de E é ponto interior a E. 

Definição 8. Em relação á W', dizemos que um conjunto E é fechado se e somente se o seu 
complementar, R" \ E é aberto. 

Definição 9. Um conjunto E C M" chama-se compacto quando o mesmo é limitado e fechado. 

Teorema 7 (Weierstrass). Seja E um conjunto compacto de M". Então, toda função contínua 
f : E ^ M assume um valor máximo e um valor mínimo. Ou seja, existem xq,xi G E tais que 
f{xo) < f{x) < f{xi) para todo x £ E. 

Prova: Ver [LimlO], p. 45. O Teorema de Weierstrass também é válido para espaços compactos 
quaisquer [MunOO]. 



Capítulo 3 

EDO's - Tópicos em Teoria Qualitativa 



Para estudar sincronização e sua consequente estabilidade, é necessário impor que os elementos 
da rede estudada possuam dinâmica. Em particular, essa dinâmica pode ser uma dinâmica não- 
linear (o modelo de rede que iremos estudar será formalmente apresentado no Capítulo 6). 

Encontrar as soluções de uma equação diferencial não-linear pode não ser uma tarefa fácil. 
Na maioria dos casos (principalmente com sistemas em altas dimensões) pode não ser possível 
de se computar analiticamente as soluções [BN69]. Na verdade, existem apenas alguns poucos 
casos em que somos capazes de encontrar analiticamente a solução [soluções] de uma equação 
diferencial [sistema de equações diferenciais], como por exemplo em equações diferenciais lineares 
cujos coeficientes são constantes. Portanto, neste capítulo trabalharemos alguns pontos da teoria 
qualitativa de equações diferencias ordinárias (EDO's). 

Seja D um aberto de M™", f : D ^ M™, m > 1, um campo vetorial de classe C", r > 2, e 
considere a equação diferencial 



Para um dado x{s) £ D, o problema de valor inicial é encontrar um intervalo I B s, da reta 
real, onde s é ponto interior e x{t) está definida em / satisfazendo a equação (3.1) e respeitando 
a condição inicial x{s) = Xg- Como / é de classe C"~, r > 2, garante-se o teorema de existência e 
unicidade das soluções (Teorema 21). Na maioria dos casos, vamos considerar / C M+ e s = 0. 

Por vezes, nos referimos â equação diferencial (3.1) pela palavra sistema e / pela palavra dinâ- 
mica. Muitos nomes são dados â solução de uma equação diferencial, fazemos a escolha de usar a 
palavra trajetória ou simplesmente solução. 

Definição 10. Seja (p : I x D ^ M*" a aplicação (j){t,XQ) = x{t). Dizemos que 4> é o fluxo do 
campo de vetores f onde x(í) : / — )• é a trajetória de f por xq em cada ponto t £ I. 

Podemos utilizar também a notação (pti^o) para o fluxo. De forma mais simplificada, o fluxo é 
uma aplicação que mapeia a condição inicial na solução. Escrevendo (/>t(x(s)) = (f)t,s{x) e conside- 
rando a operação composição, nota-se que o fluxo é um grupo de difeomorfismos a um parâmetro 
[Arn92], onde são válidas as propriedades 

1. (j)t^s ° <Ps,u{x) = (pt,uix) para todo s,u,t G I e x G D. 

2. (j)t,t{x) = X para todo t G I e x £ D. 

Exemplo 2. Dado o sistema linear x = Ax, onde A G Mat(R, m), temos que para cada ponto 



é a solução de x = Ax com a condição inicial x{tí) = xq, de modo que o fluxo çí) : M x — M™ é 
dado por 

(t){t,xo) = (f)t{xo) = exp{tA)xo. 

Para t = O o fluxo resume-se ao operador identidade e para cada t, (pt = exp(í^) é um isomorfismo 
linear cujo inverso é o fluxo (pt = exp(— ívl). 





Xq £ dado, 



x{t) = exp{tA)xQ 
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3.1 Sistemas Dissipativos 

Definição 11. Dizemos que um conjunto Q C D é invariante pelo fluxo (j) do campo de vetores f 
se 4>t{^) Ç Ç] para todo í G M. 

Equivalentemente, dizemos que íl é positivamente invariante se 0í(í2) Ç Çl para todo t > 0. 
E negativamente invariante quando (pti^) Q ^ para í < 0. Em outras palavras se x{s) £ Í7 
então x{t) G para todo í > s > O, no caso de Çl ser positivamente invariante. De forma mais 
intuitiva, dizemos que o conjunto é positivamente invariante se as trajetórias entrarem em íl lá 
permanecendo por todo o tempo futuro. Assim, temos a 

Definição 12. Dizemos que o sistema (3.1) é dissipativo se suas soluções entram, a tempo finito, 
em um conjunto positivamente invariante O, C D. 

Í7 será chamado de domínio absorvente do sistema x = f{x). A existência de um domínio 
absorvente garante que as soluções do mesmo são limitadas e portanto existam globalmente em 
decorrência do teorema de extensão (Teorema 22). A seção a seguir tratará de condições para 
garantir a existência do domínio absorvente. 



3.1.1 Segundo Método de Lyapunov 

Esta seção é de valorosa importância pois a técnica que traz a ideia da função de Lyapunov 
será usada principalmente para garantir a existência global das soluções e encontrar regiões de 
estabilidade assintótica. 

Descoberta por Aleksandr Mikhailovich Lyapunov (1857 - 1918) no final do século 19, a técnica 
é conhecida por segundo método de Lyapunov ou método direto [BN69], pois a mesma pode ser 
aplicada diretamente às equações diferenciais sem ter nenhum conhecimento das soluções. A função 
de Lyapunov é bastante utilizada para estudar a estabilidade de pontos fixos. Porém, esse não será o 
nosso propósito no uso desta ferramenta. A estabilidade de pontos fixos, será estudada no Capítulo 
4. A dificuldade desta técnica é que não existe uma fórmula fechada para encontrar tal função para 
o sistema de equações estudado. 

Seja um subconjunto conexo não- vazio de M"*. Seja V : — )• M uma função escalar conti- 
nuamente diferenciável. Definimos as noções de definitude^ de V: 

Definição 13. Dizemos que a função escalar V : M™ — )• R é positiva- definida com relação ao 
conjunto uo C M™, contendo a origem, se V{x) > O para todo x £ M™ \ uj e V{0) = 0. 

Definição 14. Dizemos que V : M™ — )• M é negativa- definida em relação ao conjunto B G 
se —V é positiva- definida com respeito à esse conjunto. 

A derivada de V com respeito ao sistema (3.1) é o produto interno 

V'{x) = {VV{x),f{x)), (3.2) 

onde 'VV{x) = (^§^{x), • • • , -^—{x)^- De fato, visto que, se x{t) é solução da equação (3.1), então 
pela regra da cadeia tem-se que 

= {VV{x),f{x)). 

Definição 15. Sejam f : D ^ M™" um campo vetorial no aberto D Ç M™ e V : — t- M uma 
função contínua. Dizemos que a função V é uma função de Lyapunov para f com respeito ao 
conjunto C D se 



^Classe dos conceitos de função positiva-definida e negativa-definida 
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(i) V é positiva- definida em relação ao conjunto e 
(ii) V' é negativa- definida em relação ao mesmo conjunto Çl. 
Definição 16. Dizemos que V é radialmente ilimitada se 

lim V{x) = oo. 

Definição 17. Dizemos que a função escalar continua V : — t- M é própria se dado um conjunto 
compacto G M então a pré-imagem V~^{K) é compacto de M*". 

Lema 2. Uma função contínua V : M™ — t- M é própria se e somente se é radialmente ilimitada. 
Prova: Ver [Tcr09], p. 180. 

Corolário 2. V é radialmente ilimitada então os conjuntos de nível V{x) = c são compactos. 

Teorema 8 (Lyapunov). Seja V uma função de Lyapunov para f em (3.1) com respeito ao conjunto 
O, C D contendo a origem, e além disso suponha que V é radialmente ilimitada. Sendo assim, as 
trajetórias da equação (3.1) entram, a tempo finito, em íl, e lá permanecem por todo o tempo futuro. 

Dizer que as trajetórias entram a tempo finito no conjunto absorvente Çl significa dizer que o 
sistema (3.1) é dissipativo. Certamente o Teorema 8 depende da existência da função escalar V 
com as propriedades citadas. Porém tal resultado não nos mostra como construir a referida função, 
e esta é a principal limitação deste método. Então, não existe um procedimento específico para a 
construção de V, mas em boa parte dos casos, uma forma quadrática é uma boa candidata. Iremos 
então considerar apenas as formas quadráticas para a função de Lyapunov e assumir que a mesma 
é dada considerando a 

Suposição 1. Existe uma matriz positiva- definida Q tal que 

V{x) = l/2{x-a,Q{x-a)) (3.3) 
onde a G M™" é fixo, V é positiva- definida em relação ao conjunto 

í): ={x G M"^ : F(a;) < p}, (3.4) 
para algum p £ e V'{x) é negativa- definida em relação à íl. 

O papel de a G M™' é realizar uma translação de eixos, se necessária. A existência de um domínio 
absorvente para o sistema (3.1), garante que as soluções são limitadas e existem globalmente, pois 
como construído, o domínio absorvente O, é um cojunto compacto de forma que o resultado segue 
pelo Teorema de Extensão (22). Vamos abordar apenas heuristicamente a prova do Teorema 8. Uma 
prova formal poderá ser encontrada em [BN69]. Vamos olhar para esse teorema geometricamente 
(veja a Figura 3.1 para o caso m = 2). De forma específica, vamos discutir a condição V'{x) < O, 
onde V é positiva-definida com respeito a uma região Í7 de M*". 

Seja c uma constante positiva e considere a equação V{x) = c. Essa equação define superfícies 
de nível em M*", as quais indicamos da forma 

Sc = {xeR"' : V{x) < c}. 

Como V é radialmente ilimitada, tais conjuntos são compactos e portanto positivamente invariantes. 
Considere que a trajetória x(í), solução de (3.1), está fora de e que Sc D íí. Por definição, 
V'{x) = {W{x), f{x)) onde o vetor 'W{x) é um vetor normal a superfície de nível V{x) = c. 
Então, a hipótese de que V'{x) < O significa que f{x) deve apontar para o interior da superfície de 
nível Sc- Por outro lado, o vetor f{x) é um vetor tangente a trajetória do sistema x = f{x) em cada 
ponto x{t). Portanto, a trajetória cruza do exterior para o interior da região limitada pela superfície 
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de nível V{x) = c para todo c. Isso faz com que quando c tenda ao bordo de Í7 as trajetórias entrem 
em Í7. Uma vez em íl, as trajetórias não deixam Q,. Suponha que x{t) G díl. para algum í = s > O 
e que x{p) ^ Í7 para p ligeiramente maior que s. Esta hipótese contradiz o fato de que V{x{t)) é 
descrescente ao longo da trajetória ^. 



VV(x) 




Figura 3.1: Ilustração geométrica para o Teorema 8 com m = 2. 



3.1.2 Sistema de Lorenz 

Um exemplo de dinâmica não-linear é a conhecida dinâmica de Lorenz. Utilizaremos tal dinâmica 
para efetuar as simulações computacionais dos principais resultados que serão apresentados ao longo 
da presente dissertação. 

Edward Lorenz (1917 - 2008) foi um meteorologista do MIT interessado em previsões do clima 
à tempos longos. Com o advento do computador, muitas pessoas direcionavam suas tentativas de 
previsão do clima a partir da análise numérica de equações que governam a dinâmica da atmosfera. 
Algo que era muito usado eram as aproximações estatísticas, especialmente a regressão linear. 
Lorenz, porém, acreditava que tais métodos eram essencialmente falhos, pois as reais equações de 
evolução estavam longe de serem lineares. Depois de experimentar vários exemplos para testar suas 
ideias, Lorenz ao estudar um trabalho por B. Saltzman [Sal62], concernente a convecção de fluídos 
térmicos, percebeu que o sistema de equações estudado por Saltzman era o ideal para provar o seu 
ponto de vista [ViaOO] . A maioria dos modelos de previsão para o clima envolve equações diferenciais 
parciais, porém Lorenz encontrou um modelo muito mais simples de ser analisado [HSD04]. 

O sistema de Lorenz, x = /(x), é tal que 

/ o-(y-x) \ 

/(x) =lx{r-z)-y\, (3.5) 
\ -bz + xy ) 

com x = {x,y,z) G M^; cr, r e b constantes positivas, conhecidas respectivamente por número de 

^Note que V'{x(u))du = V{x{t)) — V{x(s)) < O para todo í > s, portanto V{x{t)) é descrescente ao longo da 
trajetória x{t). 



3.1 



SISTEMAS DISSIPATIVOS 13 



Prandtl, número de Rayleigh, e uma proporção física, e seus valores clássicos são: cr = 10, r = 28, 
e b = 8/3. Esses são os valores que iremos considerar sempre que trabalharmos com o sistema de 
Lorenz. 

Lorenz, ao utilizar condições iniciais distintas, porém muito próximas, para este sistema, veri- 
ficou que as trajetórias divergiam gradativamente, de forma a produzir resultados bem diferentes, 
ou seja, ele percebeu que tal sistema apresenta uma dependência sensível das condições iniciais, e 
essa é a essência do comportamento caótico. Dessa forma, assumindo que o clima tem um compor- 
tamento parecido com este modelo, Lorenz concluiu ser impossível fazer previsões à longo prazo. 
Resumidamente, o sistema de Lorenz exibe uma dinâmica caótica [ViaOO]. 

Propriedades do Sistema de Lorenz 

A não-linearidade do sistema de Lorenz está associada aos termos xy e xz. Como o sistema 
de Lorenz é um sistema caótico, as trajetórias, para pequenas diferenças nas condições iniciais, 
certamente irão divergir ao longo do tempo. Considere duas equações diferenciais, ambas com 
o sistema de Lorenz, a saber, xi = /(xi) e X2 = /(x2), onde xi(í) = (xi(í), e 
X2(í) = {x2{t) , y2{t) , Z2{t)) com condições inciais próximas, porém distintas. A partir de integra- 
ções numéricas, utilizando o método Runge-Kutta de quarta ordem, podemos verificar o fenó- 
meno que aqui expomos. Utilizando as seguintes condições iniciais para o sistema: (—7, 10, 5) e 
(—7.01,10.01,5), exibimos as séries temporais xi{t) e X2{t) que representam as primeiras compo- 
nentes das trajetórias para o sistema considerado com as respectivas condições iniciais indicadas. 
O resultado é apresentado na Figura 3.2. 

20 
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Figura 3.2: O eixo horizontal representa o tempo e o vertical as componentes Xi{t) e X2{t) das trajetórias 
cuja diferença das condições iniciais, considerando a norma euclidiana, é aproximadamente 0,014. 

Temos o mesmo sistema com diferença nas condições iniciais bem próximas, entretanto, essa 
diferença cresce com o tempo porém a mesma é limitada. Mesmo considerando condições iniciais 
bem distintas, por exemplo, (0,2,0) e (O, —2,0), ambas as trajetórias divergem e se assemelham à 
borboletas confinadas em uma mesma região do (veja Figura 3.3). Isto nos elucida um importante 
fato sobre o sistema de Lorenz: todas as soluções, que não são pontos de equilíbrio, tendem ao mesmo 
conjunto. Tal conjunto é chamado de atrator estranho de Lorenz [ViaOO]. 

Função de Lyapunov Para o Sistema de Lorenz 

Vamos discorrer, nesta subseção, sobre a construção da função de Lyapunov para a dinâmica 
de Lorenz, e por consequência conseguimos também o domínio absorvente para a mesma. 
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Figura 3.3: Típica trajetória para o sistema de Lorenz considerando os parâmetros a — 10, r — 28 e 
b = 8/3. 

Considere a função 

V{x) = l/2(x-a,Q(x-a)) 

onde X = (x, y, z), a = (O, O, 2r) e Q = diag(r, a, a), a qual deve ser positiva-definida com relação a 
algum C M^. De forma explícita a função V se lê 

y{x, y, z) = l/2((x, y,z- 2r),Q{x, y,z- 2r)) 

= l/2((a:, y,z- 2r), (rx, ay, a{z - 2r))) 
= l/2[rx^ + cjy^ + a{z - 2rf]. 

Note que a condição V{x, y, z) = O O define um elipsóide em M'^ com centro em (O, O, 2r). Vamos 
então mostrar que as soluções de (3.5) entram, e permanecem confinadas, dentro da região limitada 
pelo elipsóide rx'^ + cry^ + a{z — 2r)^ = 2c para algum c > 0. Defina 

n = {{x, y, z) G : rx^ + ay^ + a{z - 2rf < 2c}. 

A derivada da função V se lê 

, dx dy dz 

V{x,y,z)=rx-+ay-+a{z-2r)- 

= raxy — rax^ + raxy — axyz — ay"^ — haz^ + 2rbaz + axyz — 2raxy 
= —rax^ — ay^ — ba{z'^ — 2rz) = —a \rx^ + + ^(^^ — 2rz)\ , 

mas h{z^ — 2rz) também pode ser escrito como b [(z — r)-^ — r^] de forma que temos 

V'{x, y, z) = -a [rx^ + + b{z - rf - br'^] . (3.6) 

A equação rx^ + y^ + b{z — r)^ = p também define um elipsóide quando p > 0. Neste caso, 
teremos V' < O quando p > 6r^. Além disso, para que sejam satisfeitas simultaneamente as duas 
condições para que V seja uma função de Lyapunov devemos escolher c > O tal que o elipsóide 
rx^ + y^ + b{z — r)^ = òr^ esteja inteiramente contido no elipsóide rx^ + uy^ + a{z — 2r)^ = 2c. 
Portanto, temos dois elipsóides 

x^ y^ (z- r)^ ^ x^ y^ (z - 2rf ^ 

^' br^ br^ r2 ^ ' (2c/r) ^ (2c/a) (2c/a) 
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Então, devemos procurar o valor mínimo de c para que se cumpra a condição Ei C -£'2- Isso 
pode ser feito através do método dos multiplicadores de Lagrange sobre dada a restrição V < 0. 
Considerando os valores clássicos dos parâmetros do Lorenz, ou seja a = 10, r = 28, b = 8/3, o 
valor de c que devemos tomar é 

Para maiores detalhes veja [Spa82]. 

Assim, podemos redefinir o conjunto Q da forma 

n = {{x, y, z) : rx^ + ay^ + g{z - 2rf < h'^r'^/{h - 1)}, (3.7) 

o qual é uma região limitada e fechada de M'^, portanto um conjunto compacto. Assim, como são 
satisfeitas as condições do Teroema de Lyapunov (8), qualquer solução do sistema de Lorenz que 
tem condições iniciais fora de Çí convergirá kVt e permanecerá dentro de por todo o tempo futuro. 
Podemos então, enunciar a seguinte 

Proposição 2. O sistema de Lorenz é um sistema dissipativo e suas trajetórias certamente entram, 
a tempo finito, no domínio absorvente (3.7). 



Capítulo 4 



Estabilidade em EDO's Lineares 
Não- Autónomas 

Vamos estudar algumas noções acerca de equações diferenciais lineares não- autónomas, em 
especial, a estabilidade de seu ponto fixo trivial. Nosso objetivo é transformar o problema de estudar 
a sincronização e sua estabilidade num problema de estudar a estabilidade da solução trivial de uma 
equação diferencial linear não-autônoma. 

Considere A{t) : J — t- Mat(M, n) uma função matricial contínua e uniformemente limitada que 
toma valores t € J M_|_ e a equação diferencial linear não-autônoma 

x = A{t)x. (4.1) 

A equação (4.1) satisfaz as condições de existência e unicidade dada a condição inicial x{s) = Xg e 
mais do que isso, as soluções estão definidas para todo t G M_|_ [Tesl2]. 

Note que x{t) = O para todo t S M+ é uma solução de (4.1) para um dado problema de valor 
inicial. Chamamos tal solução de solução trivial. Por todo este capítulo focaremos o estudo da 
estabilidade da solução trivial de (4.1) com condição inicial x{s) G M™". 

Para equações diferenciais lineares autónomas, a estabilidade de x = O pode ser caracterizada 
em termos dos autovalores da matriz constante A{t) = A. Quando, porém, consideramos que as 
equações diferenciais lineares são não-autónomas, os autovalores já não servem mais como critério de 
caracterização da estabilidade [Cop78]. De fato, podemos ter casos em que a matriz A{t) tem todos 
os autovalores com parte real negativa e ainda assim, a solução trivial seja instável, ver [Cop78] 
para exemplos. Ou ainda, casos em que as soluções não-triviais de (4.1) convergem ao vetor nulo, 
porém uma pequena pertubação pode destruir sua estabilidade, ver [BV07] e [JR08] para exemplos. 

4.1 Alguns Tipos de Estabilidade 

Para tratar os principais resultados aqui expostos, necessitamos de algumas definições. A teo- 
ria de estabilidade é ampla e vários conceitos sobre estabilidade são encontrados [Tesl2], [DL07]. 
Abordamos alguns desses conceitos que servirão ao nosso propósito. 

Definição 18. Dizemos que o ponto fixo x = O de (4.1) é estável no sentido de Lyapunov em 

t = s se para todo e > O existe õ = ô{s, e) tais que 

\\x{s)\\ <6 \\x{t)\\ < e V í > s. 

Assim, a estabilidade segundo Lyapunov é definida sobre um tempo inicial e, em particular, 
esse tipo de estabilidade não requer que as trajetórias que começam próximas da solução trivial 
convirjam assintoticamente à mesma. Outro importante conceito de estabilidade é a estabilidade 
uniforme. E um caso especial de estabilidade segundo Lyapunov, que garante que o ponto fixo não 
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esteja perdendo a estabilidade. Assim, a diferença de tal conceito com relação ao definido em 18 é 
que 6 não depende do tempo inicial s. 

Definição 19. Dizemos que o ponto fixo x = O de (4.1) é assintoticamente estável em t = s se 

o mesmo é estável no sentido de Lyapunov e se é localmente atrator, isto é, existe uma vizinhança 
U{0) 3 x{s) tal que x{t) converge assintoticamente à O quando t — )• oo, com x(t) G U{0) para cada 
t. Ou seja 

lim x{t) = 0. 

f— >oo 

Estamos, porém, interessados na estabilidade uniformemente assintótica. O interesse nesse tipo 
de estabilidade está baseada no fato de que a mesma é preservada quando o sistema sofre pequenas 
pertubações caracterizando assim a propriedade de persistência. 

Definição 20. O ponto fixo x = O de (4.1) é uniformemente assintoticamente estável se 

existe 6 > O, independente de t, tal que, para todo e > O, existe r = r(e) satisfazendo ||a;(s)|| < 6 e 
||a;(í)|| < e para todo t > s + t. 

4.2 Operador de Evolução 

Para o nosso propósito, precisamos de uma forma fechada para expressar as soluções de (4.1). 
A partir da teoria de equações diferenciais ([Tesl2]) nós podemos escrever tal solução sob a forma 

x{t) = T{t,s)x{s) (4.2) 

onde T{t,s) é chamado de operador de evolução e x{s) é a condição inicial. O operador de 
evolução satisfaz as seguintes propriedades: 

1. T{t,s)T{s,u) = T{t,u) 

2. T{t,s)T{s,t) = Im 

A propriedade 1 é imediata tendo-se em vista que podemos escrever T{t, s)x{s) = (ptix^s)) = 
4)t,s{x) e as propriedades sobre o fluxo (ver página 9). Considerando a propriedade 1, podemos 
sempre escolher t = u, produzindo T(u, s)T{s,u) = T{u,u) = Im, donde segue que T{t,s) é um 
isomorfismo com inversa dada por T~^(t,s) = T{s,t), valendo-se assim a propriedade 2. 

Exemplo 3. Quando o sistema (4.1) é autónomo, ou seja, quando A{t) = A, o operador de evolução 
do mesmo se lê 

T{t,s) =exp[(í-s)A]. 

Existem processos para se encontrar explicitamente o operador de evolução (ver [Tesl2] por 
exemplo), porém não estamos interessados em obter analiticamente tal operador mas em estimar 
sua magnitude. A seguir, daremos uma definição necessária à abordagem dos principais resultados 
aqui expostos. 

Definição 21. Seja T{t,s) o operador de evolução associado ao sistema (4.1). Dizemos que tal 
operador possui contração uniforme se 

\\T{t,s)\\ < ke-'^^*-'^ 

com k e f] positivos. 

O termo contração uniforme, associado ao operador de evolução, significa que tal operador 
de fato possui contração, como expressado pela norma, pois a mesma converge exponencialmente 
rápido à zero, e que as constantes r] e k não possuem dependência do tempo inicial s. 

O próximo teorema à ser abordado, relaciona intimamente as definições 20 e 21. 
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Teorema 9. Seja T{t,s) o operador de evolução da equação (4.1) e considere a solução trivial da 
mesma x{t) = 0. Tal solução é uniformemente assintoticamente estável se e somente se T{t,s) 
possui contração uniforme. 

Prova: Faremos primeiro a volta. Escrevamos a solução de (4.1) sob a forma x(t) = T{t, s)x{s) e 
então tomamos a norma em ambos os lados, assim 

\\x{t)\\ = \\T{t,s)x{s)\\ 

<\\Tit,s)\\\\xis)\\ 
< ke-'^^'^''>\\x{s)\\. 

Queremos mostrar que x{t) = O é uniformemente assintoticamente estável, ou seja, deseja-se que 
exita (5 > O tal que para todo e > O, exista r = r(e) > O tais que, sendo ||x(s)|| < (5 =^ 112^(^)11 < 
e Ví > s + r. Dessa forma, devemos ter 

-r]{t - s) < In (e/ôk) 

t > s + (l/r/) ln{6k/e). 

Assim, basta tomar r = (l/?/) In (ôk/e) para garantir o que desejamos. 

Para provar a ida, considere que a solução x(í) = O é uniformemente assintoticamente estável. 
Então existe (5 > O tal que para todo e > O, exite r = r(e) > O, com ||x(s)|| < 6, e ||x(í)|| < e, onde 
t > s + T. Escolha então e = (5/c, c > 1, e considere a sequência tn = s + nr, com n £ N. Perceba 
então que 

= ||r(í,.)x(.)||<^ 



para qualquer ||x(s)||/(5 < 1, o que implica 

||r(í,s)|| = sup ||r(í,s)u|| < -, 

||u||<l C 

onde u = \\x{s)\\/õ. Pela propriedade 1 do operador de evolução, ou seja, T{t,u) = T{t, s)T{s,u), 
podemos fazer 

||r(Í2,s)|| = \\TÍs + 2t,s + t)T{s + t,s)\\ 

< ||r(s + 2r,s + T)||||r(s + r,s)|| 
1 



Assim, por indução, temos 



\\ntn,s)\\<j^ 



onde, podemos impor que l/c" = e Então —rj{tn — s) = Inc mas tn — s = nr, o que 

implica que podemos tomar 77 = In {c)/t, de forma à conseguir que 

\\T{tn,s)\\<e-^^'--'l 

Fica ainda um último ponto a resolver, que é exatamente que neste caso, temos uma contração 
uniforme para o operador de evolução, porém escrita em tempo discreto. Para torná-lo em tempo 
contínuo, utilizamos o Teorema do Valor Intermediário e consideramos o caso geral em que t = 
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s + V + nr, onde O < v < t. Assim 

||T(í,s)|| < e-^í*"-^) 

= g-v{t-s-v) 

com k = e'^'". ■ 

4.3 Persistências da Solução Trivial 

Dizemos que o sistema x = f{t,x), f : M-|- x M™ — )• M™, é persistente, quando o mesmo mantém 
suas propriedades qualitativas sob pequenas pertubações ou mudanças no campo de vetores /, veja 
no livro [GH83] uma discussão mais detalhada. Mas, o que entendemos por "pequenas perturbações"? 
Muitos sistemas físicos interessantes possuem, tipicamente, uma dependência de parâmetros, os 
quais aparecem nas definições das equações como é o caso do sistema de Lorenz (3.5). Quando 
esses parâmetros são variados podem acontecer mudanças na estrutura qualitativa das soluções 
de forma que os pontos fixos podem perder sua estabilidade a medida que essa variação ocorre - 
essas mudanças são chamadas de bifurcações. Então essa leve mudança nos parâmetros é um dos 
exemplos do que pode ser uma pequena perturbação no campo de vetores. 

No caso particular de nosso interesse temos que / é um campo linear com /(í, 0) = 0. Seja 
g{t,x) uma perturbação de / tal que g{t,0) = O e considere o sistema perturbado y = f{t,y) onde 
f = f + g- Note que y = O é sempre solução do sistema perturbado. A pergunta então é: sendo 
X = O estável, y = O é também estável? Sob quais condições a estabilidade é preservada? Esta 
é uma pergunta interessante pois pode acotecer que tenhamos uma perturbação muito pequena e 
mesmo assim, a estabilidade da solução trivial seja destruída, veja [BV07] para exemplos. Diremos 
portanto que a estabilidade de x = O em (4.1) é persistente quando a mesma é preservada sob 
pertubações pequenas no campo de vetores. 

Voltando agora à teoria das contrações uniformes, sabemos que a mesma está muito relacionada 
com a propriedade de persistência da estabilidade. Os resultados desta seção tratam dessa relação. 

Teorema 10 (Persistência). Suponha que A{t) é uma função matricial contínua em um intervalo 
J C M+, referente à equação (4.1). Suponha também que o operador de evolução T{t,s) da mesma 
possui contração uniforme. Considere B{t) uma função matricial contínua em J satisfazendo 

snp\\B{t)\\=ôo < j. 

Então o operador de evolução T{t, s) da equação perturbada 

y = [Ait) + Bit)]y (4.3) 
também satisfaz ||T(í,s)|| < A;e~'^^*~'^^ com-j = rj—ôok, onder] ek são tais que \\T{t,s)\\ < ke~'^^^~'^\ 

Para a prova do teorema sobre a persistência da estabilidade da solução trivial, utilizaremos dois 
importantes resultados da teoria de equação diferenciais, a saber, o Lema de Grõnwall e o Método 
da Variação dos Parâmetros, que podem ser encontrados no Apêndice B com mais detalhes. 
Prova: Usando o método da variação dos parâmetros em (4.3), obtemos 



y{t) = T{t,s)y{s)+ í T{t,u)B{u)y{u)du. 

J s 
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Tomando a norma sobre a equação acima, e usando a desigualdade triangular, ficamos com 



\ym<\\T{t,s)y{s)\\ + 



T{t, u)B{u)y{u)du 
<ke-^^^-^y{s)\\+ f \\T{t,u)\\\\B{u)\\\\y{u)\\du 

J s 

< ke-^'^'-'^\\y{s)\\ + f ôoke-^^'~''\\y{u)\\du, 



onde as sucessivas desigualdades foram obtidas através da propriedade sub-multiplicativa da norma 
induzida (Lema 1). 

Neste ponto, tomamos a função uj{t) = e'''*||í/(í)||. Multiplicando a última desigualdade por e^*, 
obtemos ^ 

uj{t) < kui{s) + / ôç)ku}{u)du 



a qual podemos usar o lema de Grõnwall, de forma que obtemos 

u}{t) < kuj{s)e^oHt-s)^ 

Substituindo nesta última desigualdade a função uj{t), segue que 

||y(í)||<te(-''+'5o'=)(*-^)||y(5)||. 

Logo, podemos tomar ||T(í, s)\\ < ke~"'^^~^\ com J = rj — 6ok. ■ 

O próximo resultado é importante pois o mesmo dá as hipóteses sobre a qual o resto de Taylor, 
advindo da linearização pela expansão em série de Taylor, não destrói a estabilidade da solução 
trivial. 

Proposição 3 (Princípio da Linearização). Considere que o operador de evolução da equação (4.1) 
tem contração uniforme. Considere a equação perturbada 

y = A{t)y + Rit,y) (4.4) 

onde R{t, y) é tal que R{t, 0) = O, V í > O, e a seguinte propriedade é satisfeita 

V e > 0,3 5 > O : se ||y|| <õ^ \\R{t,y)\\ < e||2/||,V t > 0. 

Então o operador de evolução de (4.4) tem contração uniforme. 

Prova: Utilizando a variação dos parâmetros na equação (4.4) tem-se 

y{t) = T{t, s)y{s) + T T{t, u)R{u, y{u))du (4.5) 



s 



onde T{t,s) é o operador de evolução de (4.1), que por hipótese satisfaz ||T(í,s)|| < ke ^\ 
Calculando a norma em (4.5) e utilizando as hipóteses dadas tem-se 

\\yit)\\ < \\Tit,s)\\\\y{s)\\+ f \\T{t,u)\\\\Rin,yiu)\\du (4.6) 



< A:e-''(*-^)||y(s)|| + / efee^^^*-") ||y(n) (4.7) 
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Considere a função real uj{t) = e^*||y(í)||. Multiplicando e''* em (4.7) teremos 

Lo{t) < kuj{s) + ke oj{u)du. (4.8) 

A partir daqui, utilizamos o Lema de Grõnwall em (4.8) e seguimos os mesmos passos utlizados na 
prova do Teorema 10 mostrando que 

||y(í)|| <fce(-''+^^)(*-^)||2/(s)||. 

Escolhendo e < r///c, encerra-se a demonstração. ■ 

Algo interessante que a Proposição 3 nos diz é que a estabilidade da solução trivial de (4.1) não 
muda mesmo considerando a pertubação não-linear R{t, y) com a propriedade descrita, isso porque 
podemos tomar e tão pequeno quando desejarmos. Note ainda que a condição dada sobre i?(í, y) é 
a mesma condição dada sobre o resto de Taylor [LimlO], isto é, limy_j>o R{t,y) /\\y\\ = 0. 



4.4 Critérios para Contração Uniforme 

O resultado a seguir fornece um importante critério a respeito de como garantir que o operador 
T{t,s) de (4.1) tenha contração uniforme. O Teorema da Diagonal Dominante garante que os 
coeficientes da diagonal da matriz A{t) de (4.1) controlam a estabilidade do sistema. 

Teorema 11 (da Diagonal Dominante). Seja A{t) = [Aíj(í)]™^^ uma função matricial contínua 
e limitada com t £ M+. Suponha que existe uma constante r] > O tal que 

m 

A^^{t)+ l^y(í)l<-^<0 (4.9) 

para todo t £ M_(_ e i = 1, . . . ,m. Então o operador de evolução de (4.1) tem contração uniforme. 

Prova: Seja x = (xi, • • • ,Xm) G I^*"- Usaremos a norma ||x||oo — max^ p^-ra estimar x(í). Kixe 
u > O, então para algum 1 < i < m tem-se que ||x(íx)||oo = |2;í(u)| ou ainda ||2;('u)||^ = Xi'^{u). 
Analisaremos a norma ||x(í)||^ pois tal função é diferenciável. Note porém que poderia acontecer 
de que exatamente neste u fixado a norma ||2;(ii)||^ ser atingida em mais de um i. Suponha por 
exemplo, sem perder a generalidade, que esse máximo é atingido em h e k, h,k G {1, • • • , m}. Então 
surge a pergunta: para t numa vizinhança de u, qual seria maior, Xh'{t) ou Xfc^(í)? A resposta é: 
independente de quem seja maior, esse máximo continuará sendo o máximo por pelo menos um 
intervalo (u, u + e) = / com e > O suficientemente pequeno. Suponha que o que acabamos de 
afirmar não seja verdade, isto é, suponha que para t £ I, o máximo alterna-se entre as componentes 
h e k um número arbitrário de vezes. Neste caso, teríamos que lime^o l|3^(^ + c)lloo existe e 
portanto (i/(ií||x(í)||^ não existiria t = u, o que é absurdo pois todas as componentes Xi^{t) são 
diferenciáveis visto que x(í), solução não-trivial de (4.1), é diferenciável. Portanto podemos impor, 
sem perder a generalidade, que se Xfc^(n) = ||x(n)||^ então Xfc^(í) = ||a^(í)lloo para todo t £ I. 
Então, tomando t G I tem-se 

^|l|x(t)||L = x.(í)|x.(t). 

Mas, note que de (4.1) tem-se 

d ™ 
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assim 



Akk{t)xk^{t) + ^Akj{t)xj{t)xk{t) 



< 



Akkit) + Y\Akjit)\ 

jy^k 



Xk^{t), 



onde por hipótese vale (4.9). Então 



^"x(í)||L<-2W(í) = -2r?||xr"2 



dt 

Tomando a integral, definida em [s, t] C /, temos 

\\x{t)\\1<\\x{s)\\1-2r^ f\\x{v)\\ldv, 

J s 

da qual podemos usar a desigualdade de Grõnwall, obtendo 

||x(í)||oo<e-^(*-^)||x(5)|U 
, é estritamente negativa, então ||x 



Como a derivada de ||x 
em í E I e portanto vale que \\x 



(4.10) 

não possui máximo local 



, < ||x(s)||^ para t > s. Mais do que isso, neste caso, é possível 
garantir que ||x(í)||oo é uma função monotonicamente decrescente em M_|_ (veja [Cop78], a partir 
da página 55 para mais detalhes), então a desigualdade (4.10) será válida para todo t G M+, t > s. 
Como em Mat(M, m) todas as normas são equivalentes, então existem c > O e /c > O tais que 
c||T(í,s)||oo < ||T(í,s)|| < A:||r(í,s)||oo, portanto ||r(í,s)|| < ke-^^^-^l ■ 

Agora, enunciamos um novo critério que está relacionado com o espectro da parte simétrica da 
matriz A{t) de (4.1). 

Teorema 12 (da Parte Simétrica). Considere o sistema (4.1). Suponha que A{t) é uniforme- 
mente limitada em t G M4. e que sua parte simétrica é negativa- definida. Então o operador de 
evolução de (4.1) satisfaz ||r(í,s)||2 < e^^''~^\ onde r/ < O é uma cota superior para os autovalores 
da parte simétrica de A{t). 

Prova Considere x G e ||x||2 = {x,x). Considere a solução não-trivial x{t) de (4.1) definida em 
M+. A função ||a;(í)||2 é contínua e diferenciável. Então 



1 d_ 
2dt 



\x{t)\\l = l/2{A{t)x,x) + 1/2 {x,A{t)x) 



onde B(t) é a parte simétrica de A{t). Neste momento consideraremos a análise para cada t fixado. 

Como B é uma matriz simétrica então, pelo Teorema Espectral, seus autovetores ortonormali- 
zados formam uma base de M*^. Podemos representar x nessa base, x = YlT^i forma que 

faz-se valer 



(x, Bx) < ( CiVi, ^ CiXiVi ) 
\i=i i=i I 

< \m{x,x), 

onde Al < • • • < Xm < O são os autovalores de B. Como B{t) é uniformemente limitada em t e 
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negativa-definida podemos garantir que existe < O tal que Xm{i) < V para todo todo t G M_|_. 
Dessa forma, temos a relação 

j^\\xml<H\m\l 

à qual podemos tomar a integral, definida em [s,í] C M_|_, s < í, e utilizar a desigualdade de 
Grõnwall, produzindo 

||x(í)||2<e''(*-^)||a;(s)||2. 



Capítulo 5 

Tópicos em Teoria dos Grafos 



Abordaremos os fundamentos de teoria dos grafos para estruturar e modelar redes de elementos 
interagentes afim de fundamentar matematicamente o objetivo de estudar a sincronização e suas 
estabilidades em tais redes. Por todo o presente capítulo, focaremos o estudo na teoria espectral 
dos grafos, ou seja, o estudos dos autovalores da uma das matrizes que caracteriza a estrutura de 
conexão, a saber, o Laplaciano. Inicialmente daremos algumas definições necessárias ao longo do 
presente capítulo. Consideraremos apenas grafos finitos. 

Definição 22. Um grafo G = {V{G), E{G)) consiste de um conjunto finito de vértices V{G) e 
um conjunto, consequentemente finito, de arestas E{G), onde cada aresta consiste de um par de 
vértices. 

Quando ao referirmos que uma aresta {u, v} G E{G) vamos utilizar a notação « ~ u e dizer que 
u é vizinho de v. 

Um grafo direcionado consiste dc um conjunto dc vértices V{G) e um conjunto de arestas E(G), 
onde cada aresta consiste dc um par dc vértices ordenados. 

Definição 23. Dizemos que um grafo é não-direcionado se os pares de vértices {u, v} G E{G) 
são não-ordenados. 

Assim, em grafos não-direcionados, se u é vizinho de v então reciprocamente v é vizinho de u. 

Definição 24. Diremos que um grafo é simples quando o mesmo não possui auto-conexões ou 
múltiplas arestas entre dois vértices. 

Definição 25. O grau de um vértice i é numericamente igual ao número de vizinhos que o mesmo 
possui. Indicaremos o grau do vértice i por gi. 

Um caminho em um grafo é uma sequência de vértices distintos e conectados. 

Definição 26. Um grafo é conexo se para quaisquer dois vértices u e v existe pelo menos um 
caminho de u até v. 

De forma semelhante, dizemos que G é não-conexo se em G existe algum vértice que não pode 
ser alcançado por qualquer outro. Para todo grafo não-conexo G, sempre podemos decompô-lo da 
forma simbólica G = Gi U G2 U • • • U , se G tem r componentes conexas, de forma que os subgrafos 
Gi, G2, • • • ,Gr são chamados componentes conexas de G. 

Trabalharemos apenas com grafos que são simples, não-direcionados e conexos. Um grafo pode 
ser simples e ainda assim representar matematicamente uma rede complexa. Para que o leitor não 
faça confusão desses conceitos, note que a rede é simples como traz a definição acima, porém pode 
ter uma estrutura não regular de conexão entre seus elementos. Um grafo simples e regular é aquele 
o qual podemos identificar precisamente a regra que gera o conjunto de vértices. 

Para quaisquer dois pares de vértices u,v & V(G), a distância entre u e v é o menor caminho 
dentre todos os possíveis entre u e v. 

Definição 27. Em um grafo G, o diâmetro de G é a máxima distância sobre todos os pares de 
vértices em G. Utilizaremos a letra "d" para indicar o diâmetro de G. 
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5.1 Matrizes de Adjascência e Laplaciana 

Seja G um grafo simples e não-direcionado com n vértices. A matriz de adjacência de G, 
denotada por Aq ou simplesmente A = [^ij]"j=i, quando a indicação do tipo de grafo não for 
necessária, é uma matriz n x n definida da seguinte forma 



Aij — 



1 se i e j são vizinhos; 
O caso contrário. 



Observe que por construção, uma matriz de adjacência é uma matriz simétrica. Em termos dessa 
matriz, o grau Qí, do vértice i de G pode ser escrito como 



n 



9i — 5^ -^ij- 
j=í 

A matriz laplaciana de G, denotada por Lq ou simplesmente L = [Lij]fj^i, é uma matriz 
n X n definida como segue: 

{9i sei = j; 

— 1 se i e j são vizinhos; 
O caso contrário. 

Seja Dg a matriz diagonal formada pelos graus dos vértices de G. Então podemos reescrever a 
matriz laplaciana na forma 

Lg = Dg- Ag. 

De forma mais precisa, note que os elementos da matriz laplaciana e da matriz de adjacência estão 
relacionados da seguinte maneira 

Lij = õijgi- Aij (5.1) 

onde o Sij é o delta de Kronecker. 

A matriz laplaciana, por construção, também é uma matriz simétrica. Na seção a seguir, dis- 
corremos algumas propriedades relacionadas à essa importante matriz. 

5.2 Propriedades Espectrais do Laplaciano 

Resumimos as principais propriedades relacionadas ao espectro do laplaciano no 

Teorema 13. Seja G um grafo simples e não-direcionado e L sua respectiva matriz laplaciana. 

Então 

(a) Todos os autovalores de L são reais. 

(b) O (zero) é autovalor de L associado ao autovetor c = (c, c, • • • , c), c G M \ {0} o qual é único. 
Em particular, podemos considerar o autovetor 1 = (1, 1, • • • ,1). 

(c) L é positiva semi-definida e seus autovalores podem ser ordenados, de acordo com suas mul- 
tiplicidades algébricas, da forma 

O = Al < A2 < • • • < A„. 

(d) A multiplicidade do autovalor O é igual ao número de componentes conexas de G. 
Prova: 

(a) Sabemos que a matriz L é simétrica, ou seja, L = L^ . Considere A um autovalor dc L e v o 
autovetor associado à A. Assim, Lv = Xv, donde segue que v^Lv = Xv^v. Note que (v^^Lv)'^ = v^Lv. 
Logo, v^Lv = X^v^v = Xv^v, ou seja, A^ = A. 
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(b) Note que 



Lc = j ^ LijC, • • • , XI ^^j'^ ) = (O, • • • , 0) = Oc. 



Por outro lado, se Lv = O então 

n n n 

X LijVj = 0, X L2jí;j =0, • • • , J2 ^njVj = O 

j=i j=i j=i 

mas 

n n n n 

j Í=2 j=2 j=2 j=l 

assim, seguindo o mesmo argumento temos que 



3 

O que produz 



Y ~ '^l) + '^2i(^^i -V2) + --- + Y ^"i(^Í -^ri)=YY ^iò(Vj -Vi)=^ 



3 3 3 ^ 3 

implicando Vj = Vi Vi, j = 1, • • • , n. 

(c) Mostraremos que L é positiva semi-definida. Seja x = {xi, ■ ■ ■ ,Xn) € K". Então, queremos 
mostrar que x^Lx > 0. Note que 

JijXj 



c^Lx = Yj Yj ^ 

i 3 

= Xi {^ij ~ — Y^. ^ij^j ~ Y^. Y^. -^ij-^j 



13 13 ^ 3 



mas 



Então 



^Yj "^^ ^Yj ^^3-'^3 ~ ^Yj ~ ^Yj ^Yj ^í.?"^' 

i j i i j 

^ Lx = Y^, Y^. AijXi — Y^, Y^. AijXiXj, 



X 

i 3 13 

mas como A é simétrica tem-se que 

Y^j Y^j ^ij^i = Yj Yj ^3i^3 ~ Yj Yj ^^3^3 ' 
i 3 i 3 i 3 

assim, 

2x'^ Lx = Y^ Y^ AijXi^ + Y^ Y^ ^ij^j^ " ^ Yj Yj ^iJ^i^j 
i j i 3 i 3 

= ^Y^ Yj -^íÁ^í^ ~ 2xiXj + Xj"^) = Yj ^Yj ~ ^i)^ — ^ 

i j i 3 

portanto, 

x^Lx > 0. 

Como L é positiva semi-definida, segue que seus autovalores são não-negativos e portanto podemos 
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ordená-los da forma 

O = Al < A2 < • • • < A„. 

Para provar a afirmação (d), devemos notar que se o grafo G tem r componentes conexas, 
sendo Gi, . . . ,Gr suas componentes, então é possível representar L como uma divisão em blocos 
Li, . . . ,Lr. Seja m a multiplicidade algébrica do 0. Então cada Li tem um autovetor Cj associado 
ao autovalor 0. Note que Cj = (cji, Cj2 • • • , Cjn ) pode ser definida como 

J 1 se j pertence a componente i 
1^ O caso contrário , 

portanto m > r. Como o autovetores c, são únicos então m = r. M 

Analisar o espectro do Laplaciano é de grande importância ao cumprimento do objetivo de 
estudar a estabilidade da sincronização em redes. De forma mais precisa, sabemos que o segundo 
autovalor do laplaciano está intimamente relacionado ao valor do parâmetro crítico de acoplamento 
entre os elementos da rede para a sincronização da mesma e sua respectiva estabilidade. 



5.3 Espectro de Alguns Grafos Regulares 

Vamos agora examinar os autovalores e os autovetores do Laplaciano de grafos com estrutura 
regular de conexão. De forma mais precisa, vamos examinar 

1. O grafo complete de n vértices, Kn, tal que E{Kn) = {{u,v) : u ^ v}; 

2. O grafo star de n vértices, Sn, tal que E[Sn) = {(l^^^) : 2 < f < n}; 

3. O grafo path de n vértices, P„, tal que E{Pn) = {{u, ii-|-l):l<'u<n}e 

4. O grafo ring de n vértices, i?^, que além de ter todas as arestas do grafo path possui também 
a aresta (1, n). 

Proposição 4. O Laplaciano de Kn tem autovalor O (zero) com multiplicidade de 1 e autovalor n 
com multiplicidade n — 1. 

Prova: Como o grafo Kn é conexo, segue pelo Teorema 13 (d) que o autovalor O (zero) tem 
multiplicidade 1. Seja v = {v{l),--- ,v{n)) um vetor qualquer, não-nulo, ortogonal ao vetor 1 = 
(!,••• , 1), então 

i 

Assuma sem perder a generalidade que v{l) 7^ 0. Vamos calcular a primeira coordenada de Lk„v. 
Ficamos com 

n 

{LK„v){l)=giv{l)-Y,Hi) 

1=2 

n 

= (n — l)f(l) — ^^1^(0 = nv{\)^ 

i=2 

onde gi = -- -=gn = n— léo grau do vértice 1. Estendendo a análise para todas as coordenadas 
de V temos que Lk„v = nv. M 

Lema 3. Seja G um grafo conexo de n vértices e sejam i,j G {2, • • • ,n} vértices de grau 1 tais que 
ambos estão conectados à um outro vértice k. Então o vetor v{u) dado por 

{1 para u = i; 

— 1 para u = j; 

O caso contrário, 
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é um autovetor do Laplaciano de G com autovalor 1. 



Proposição 5. O grafo Sn tem autovalor O (zero) com multiplicidade 1, autovalor 1 com multipli- 
cidade n — 2 e autovalor n com multiplicidade 1. 

Prova: Como 5„ é um grafo conexo então pelo Teorema 13 {d) segue que o autovalor O tem 
multiplicidade 1. Aplicando o Lema 3 aos vértices i ei + 1 encontramos n — 2 autovetores linearmente 
independentes. Afim de encontrar o último autovalor utilizamos o fato de que o traço de uma matriz 
é igual à soma de seus autovalores. Notamos que ti{Ls„) = 2n — 2. Como identificamos n — 1 
autovalores cuja soma é n — 2, então o último autovalor é igual a n. ■ 

Antes de enunciar a próxima proposição vamos discorrer algumas palavras sobre o Laplaciano 
de Rn- De forma explícita o mesmo se lê 



2 


-1 










-1 


-1 


2 


-1 













-1 


2 


-1 















-1 


2 


-1 


-1 










-1 


2 



Dado um vetor de f = (f i, ^2, 



,Vn) S K", notamos que a ação de L/j,^ sobre v produz 
(Lr^v)^ = -Vk+i - Vk-i + 2vk, 



para \ < k < n. Note ainda que dada uma função real f{t) de classe C^, a sua segunda derivada 
pode ser escrita da forma 



lim 



f{t + h) + f{t-h)-2f{t) 

/l2 



Sendo assim, vemos que o resultado da ação de Lr^ sobre í; se assemelha a segunda derivada discreta, 
escolhendo h = 1. As funções reais tais que suas segunda derivada é um múltiplo de si mesmo são 
os senos, cossenos e exponenciais, que por sua vez pode ser expressa em função de senos e cossenos. 
Podemos então enunciar a seguinte 



Proposição 6. O grafo Rn tem autovetores 

Xk{u) 

Vkiu) - 



sm 



cos 



/ 2kUTT \ 



com o < k < n/2 se n é par e O < k < (n — l)/2 se n é ímpar e 1 < u < n. Ambos os autovetores 
têm autovalor 2 — 2cos{2kiT/n). O vetor xo{u) = O deve ser desconsiderado e se n é par então 
da mesma forma o vetor Xn/2Íu) = O deve ser desconsiderado. Note que yo{u) = 1 é o autovetor 
constante. 

Prova: Para encontrar o autovalor vamos considerar apenas o vértice 1, ou seja, u = 1 sem perder 
a generalidade. Então, sabendo que sin(2Qí) = 2sin(a) cos(a) e que cos(2qí) = 2cos'^{a) — 1, temos 

Ln„xk{l) = 2xk{l) - Xfe(O) - 2xk{2) 

= 2 sin {2k7r/n) — 2 sin (2/c7r/n) cos {2k'K /n) 

'2kTT" 



2 cos 



n 



Xfc(l). 
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Da mesma forma 



2yfc(l)-2/fc(0)-2yfc(2) 

2 cos (2fe7r/n) — 1 — cos (2 • 2k'K/n) 

2 cos {2k'K /n) — 1 — cos^ (2k'n /n) — sin^ {2k'K jn) 



2 cos 



2/í;7r 



n 



cos {2kii/n) 



2 cos 



2A;7r 



n 



Antes de enunciar o próximo resultado vamos verificar uma importante relação entre os grafos 
Pn e R2n- Vamos olhar para P„ como um quociente de i?2n pela seguinte relação de equivalência: 
Dados u Q V vértices de R2n temos que u é equivalente a u se e somente se it + f = 2n + 1. Dessa 
forma, por exemplo considerando n = 5, ou seja, considerando os grafos P5 e -Rio, em iíio estamos 
identificando os vértices 1 e 10, 2 e 9 e assim por diante, tornando o quociente de R2n pela relação 
de equivalência dada igual a P„. 

Argumentos análogos aos utilizados na Proposição 6 nos mostram que 

Proposição 7. O Laplaciano de Pn tem os mesmos autovalores de iÍ2n. e autovetores 

Vk{u) = sm{ku7r /2n + 7r/2n), 

com O < k < n. 

Teorema 14. Seja G um grafo simples de n vértices, d o seu diâmetro, gi o menor grau entre 
todos os vértices e X2 o segundo autovalor do laplaciano. Então 

4 

nd 
ngi 
n — 1 



1. A2 > 

2. A2 < 



As provas para estas desigualdades podem ser encontradas em [Moh91] e [Fie73] respectivamente. 
Agrupamos os principais resultados da análise espectral do Laplaciano na Tabela 5.1, onde gi e g^ 
representam respectivamente o menor e o maior grau no grafo. 



Grafo 


A2 


Qn 


91 


d 




n 


n — 1 


n — 1 


1 


Rn 


2-2cos(^) 


2 


2 


(n + l)/2 se n é ímpar 
n/2 se n é par 


Sn 


2 se n = 2 
1 se n > 2 


n — 1 


1 


2 


Pn 


2-2cos(^) 


1 


1 


n — 1 



Tabela 5.1: Propriedades espectrais do Laplaciano de algumas redes regulares. 



5.4 Redes Complexas 

Nesta seção abordaremos alguns exemplos de redes complexas, como elas são formandas, suas 
propriedades e aplicações. Inicialmente daremos a definição formal desse termo. 

Definição 28. Uma rede, modelada a partir de um grafo G, é dita ser complexa se G não possui 
uma estrutura regular de conectividade. 

Os exemplos que queremos abordar de redes complexas estão separados em subseções como 
segue. 
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5.4.1 Redes Aleatórias 

As redes aleatórias, ou redes que seguem o modelo Erdõs-Rényi, tem a seguinte característica 
de geração: para um dado p fixo, O < p < 1, cada aresta em potencial é escolhida com probabilidade 
p, independentemente das outras arestas. Dessa forma, todos os vértices de um grafo aleatório tem 
a mesma esperança para o grau. Nesse modelo de rede, podem existir vértices isolados, ou seja, a 
rede pode não ser efetivamente conexa. Quando o número de vértices n — )• oo então a probabilidade 
de que a rede seja conexa tende a 1 [CLotMSUOG]. 

5.4.2 Redes Pequeno Mundo 

O termo fenómeno de pequeno mundo é usado para se referir à duas propriedades distintas, 
a saber, a propriedade de pequena distância (dois estranhos estão tipicamente ligados por uma 
pequena cadeia de conhecidos mútuos) e a propriedade de efeito de agrupamento (duas pessoas 
que compartilham o mesmo vizinho possuem uma maior probabilidade de também serem vizinhos). 
Esse tipo de rede foi introduzido por Watts e Strogatz em [WS98]. Alguns exemplos desse tipo 
de rede complexa são as redes neurais, redes elétricas, redes de co-atuação de atores em filmes 
americanos, dentre outras. As redes de pequeno mundo encontram-se em um meio termo entre as 
redes regulares (não complexas) e redes aleatórias. Portanto, uma forma de gerar esse tipo de rede 
é considerando inicialmente uma rede regular tipo ring com n vértices e cada vértice com grau 
em seguida reconecta-se cada aresta de forma aleatória com probabilidade p, onde p = O gera a 
rede regular descrita e p = 1 uma rede aleatória. Essas redes apresentam, significante aumento na 
velocidade de propagação de um sinal e consequente sincronizabilidade. 

5.4.3 Redes Scale-Free 

Nem as redes aleatórias ou as de pequeno mundo tem uma propriedade frequentemente observada 
em redes do mundo real, a saber, o comportamento de que os vértices possuem probabilidade de 
ter grau g seguindo a lei de potência 

P{g) « g-^. 

Como sabemos, o grau de um vértice é o número de vizinhos que o mesmo possui. A lei de 
potência assegura que o número de vértices com grau g é proporcional a g~^ com f3 > 1. As redes 
que seguem essa lei de distribuição de potência são conhecidas como scale-free networks [DchlO]. 
Este modelo de rede é devido a Barabasi e Albert [BA99]. O termo scale está relacionado tanto 
ao espaço quanto ao tempo e de fato, escalas de espaço e tempo podem coexistir simultaneamente. 
Por exemplo, as redes de chamadas telefónicas tem formas muito similares mesmo em diferentes 
regiões geográficas e em diferentes horários [CLotMSUOG]. As redes que seguem o modelo Barabasi- 
Albert possuem uma grande quantidade de heterogeneidade, ou seja, enquanto a maioria dos vértices 
possuem apenas alguns vizinhos, alguns poucos vértices, chamados de hubs possuem muitos vizinhos. 
Um exemplo claro desse tipo de rede é a rede de estradas que conectam as cidades de um país. 
As grandes cidades possuem muitas alternativas de estradas que estão ligadas à ela, e as pequenas 
cidades poucas alternativas de ligação com outras cidades, porém o número de cidades pequenas 
é muito maior do que o de mega-cidades. Por ter essa característica de heterogeneidade, esse tipo 
de rede complexa, pode não ter sincronização global, isto é, neste caso apenas os hubs podem 
apresentar sincronização porém os demais elementos da rede ficarem fora desse estado [PcrlO]. 



Capítulo 6 



Redes de Osciladores com Acoplamento 
Difusivo 



o comportamento oscilatório desempenha um papel importante na natureza. Toda forma de 
vida exibe tal comportamento em cada nível de organização biológica com períodos que podem 
variar desde milissegundos, como é o caso dos neurónios, a anos, como é o caso da inter ação presa- 
predador na ecologia [APN99]. Considere um grafo G que modela uma certa rede. Dizemos então 
que os elementos da rede (vértices de G) possuem um comportamento oscilatório, ou chamamos 
os elementos de osciladores, quando os mesmos possuem alguma dinâmica. Um oscilador pode 
possuir qualquer tipo de dinâmica, isto é, pode ter um comportamento periódico, não-periódico ou 
caótico. No nosso caso, para qualquer que seja o tipo de comportamento que o oscilador possuir, 
representamos esse comportamento através de equações diferenciais ordinárias. 

A presente dissertação trata apenas o caso de comportamento oscilatório quando as equações di- 
ferenciais ordinárias são todas idênticas, em outras palavras, do ponto de vista isolado, os osciladores 
são todos idênticos, a menos das condições iniciais. 

Neste momento, podemos introduzir o nosso modelo de rede. Considere que um grafo G modela 
uma rede qualquer. Por padrão, consideramos sempre que n é o número de vértices em G. Então, 
para que possamos introduzir o comportamento oscilatório, sobre cada vértice da rede introduzimos 
uma cópia da Equação (3.1). A regra de nomeação dos vértices será a seguinte: considerando que 
os graus dos vértices de G podem ser dispostos da forma 



então o vértice i será aquele que possui grau Qí. O modelo de interação entre os elementos de uma 
rede, que utilizamos, é o modelo de acoplamento chamado difusivo: 

Definição 29. Dizemos que uma rede possui modelo de acoplamento difusivo quando a dinâmica 
do vértices i é influenciada por seu vizinho j de forma proporcional à uma função da diferença de 
seus estados. 

Chamaremos essa função acoplamento de H : M™ — t- IR*" e sem perder a generalidade vamos 
considerar que a mesma é uma matriz em Mat(M, m). A função de acoplamento poderia ter até 
mesmo um comportamento não-linear, porém por toda a dissertação, vamos considerar apenas o 
caso em que H é uma matriz positiva-definida. Se supormos que os vértices i e j são vizinhos, então 
a influência que o vértice j exerce sobre seu vizinho i é igual 



de modo que a análise é recíproca. Em muitos exemplos vamos considerar que H = Im, nesse caso 
esse modelo de acoplamento é chamado de acoplamento totalmente difusivo. 

Sendo assim, considerando uma rede de n osciladores idênticos, com acoplamento difusivo, tem- 
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se que a dinâmica de um vértice i qualquer da rede passa a ser descrita por 

n 

ii = f{xi) + a^Aij[H{xj) - H{x,)] (6.2) 
i=i 

onde a é o parâmetro global de acoplamento, ou seja, para todo vértice i = 1, • • • ,n temos que a 
é o mesmo valor. Note que a equação (6.2) está em função dos elementos da matriz de adjacência 
do grafo e note ainda que o modelo está bem definido pois o mesmo considera que a dinâmica do 
vértice i é influenciada apenas por seus vizinhos e de fato, Aij = O se j não é vizinho de i. 

Nos podemos também representar o nosso modelo através da matriz laplaciana do grafo, pois se 
considerarmos apenas o termo de acoplamento e lembrando que Lij = ôijQi — Aij (Equação (5.1)) 
tem-se 

Y,MH{x,) - H{xi)] = Y,A^,H{x,) - g,H{xi) 
j j 

= - ^LijH{Xj). 

j 

Sendo assim, reescrevemos a dinâmica do vértice i da forma 

n 

Xi = f{xi) - a'^LijH{xj). (6.3) 
i=i 

O modelo de acoplamento difusivo desempenha um papel importante para a sincronização de- 
vido exatamente â sua natureza difusiva, ou seja, a tentativa de igualar os estados dos vértices da 
rede, influenciando assim positivamente a rede â um estado síncrono. Observe que se todos os ele- 
mentos da rede possuem a mesma condição inicial, então os termos de acoplamento irão desaparecer 
identicamente. Além disso, se a rede em algum momento íq está em um estado de sincronização, 
então é possível mostrar que a variedade de sincronzação, ou seja 

Xi{t) = X2{t) = ■ ■ ■ = Xn{t) = S{t) (6.4) 

para todo t > to é uma variedade invariante no tempo para qualquer que seja o parâmetro de 
acoplamento a > O e para qualquer que seja a escolha da função de acoplamento H. Queremos 
caracterizar as sincronizações que são persistentes e portanto introduzimos formalmente o termo 
que é essencial nesta dissertação. 

Definição 30 (Sincronização). Dizemos que uma rede está sincronizada se existe 6 > O tal que 
para todo e > O existe r = r(e) > O satisfazendo \\xi{u) — Xj{u)\\ < 5, para algum u > O, e 
\\xi{t) — Xj{t)\\ < e para todo t > u + t, e para quaisquer i e j. 

Esta noção de sincronização produz uma variedade de sincronização que é uniformemente as- 
sintoticamente estável, o que nos garante a persistência da sincronização sob perturbações. 

O termo sincronização global, no nosso contexto, está relacionado â dois aspectos, a saber, glo- 
bal no sentido de que a sincronização não é local (o que pode acontecer, por exemplo em rede 
complexas tipo scale-free), e global no sentido de que a coincidência dos estados (Equação (6.4)) é 
preservada com o tempo t. Formaliza-se essa propriedade no Teorema 15. De modo mais preciso, 
vamos mostrar a invariância do movimento globalmente sincronizado (6.4), para isso precisamos 
reescrever o nosso modelo (6.3) em uma forma compacta onde possamos agrupar as equações do 
movimento de todos os vértices, de forma que passemos a visualizar apenas uma equação de mo- 
vimento em M"™ e não n equações em M™". Considere então X = {xi,X2,--- ,Xn) S M""™" onde 
estamos considerando o empilhamento dos vetores xi,X2, - ■ ■ ,Xn G M™. Da mesma forma considere 
F{X) = (/(xi), /(x2), • • • , f{xn))- Então o modelo (6.3) pode ser reescrito da forma 
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X = F{X)-a{L^H)X (6.5) 

onde <Si representa o produto de Kronecker. 
Sejam $í(-) é o fluxo da equação (6.5) e 

AT = {1 ^ s(í) G M"^™ : s{t) = Xi{t) G M"', V 1 < i < n} (6.6) 

a variedade em M"^™ que representa o estado síncrono global. Enunciamos então o 

Teorema 15. N é uma variedade invariante pelo fluxo. 

Prova: Vamos considerar apenas o caso em que N é positivamente invariante pelo fluxo. Dada a 
condição inicial Xq G N, isto é, Xq = 1 <^ -5(0), afirmamos que X(t) = 1 (8) s(í), í > O, é solução de 
(6.5). De fato pois 

X{t) = F{X{t)) - a{L H)X{t) 

= F{1 s(í)) - a{L (g>H){l(^ s{t)) 

= 1 /(s(í)) - a(Ll) ® Hsit) = 1 /(s(í)) 

visto que LI = O, F{1 s{t)) = 10 f{s{t)) e X{t) = 1 (8) s{t). Além disso, se Xi{t) = s{t) para todo 
i, de (6.2) segue que é = f{s). Dessa forma, o fiuxo pode ser escrito como Í>í(1 (8" •s(O)) = X{t) = 
I0s{t), logo $í(Xo) e N.m 

Esse teorema nos fornece uma forte garantia sobre o modelo que estamos utilizando, ou seja, 
ao considerar que, isoladamente, os osciladores que possuem a mesma dinâmica. Se os osciladores 
tem as mesmas condições iniciais, então os mesmos estarão automaticamente sincronizados e assim 
permanecerão por todo o tempo futuro. Se por outro lado, os osciladores não possuem as mesmas 
condições iniciais, mas para algum í > O, o parâmetro global de acoplamento garante que os 
osciladores irão sincronizar, então da mesma forma, sincronizados eles estarão por todo o tempo 
futuro. Para que esse cenário seja coerente, estamos considerando que não existem influências de 
pertubações sobre a variedade de sincronização. Porém, na prática, essa coincidência de estados só 
estará garantida se tal regime é estável, ou seja, se esse subespaço permanece inalterado após uma 
pequena pertubação [RP03]. 

6.1 Parâmetro Global de Acoplamento 

Analisaremos agora o quão importante é o parâmetro de acoplamento a para a sincronização e 
sua estabilidade. 

Considere que o campo de vetores / da dinâmica isolada dos osciladores está associado â uma 
dinâmica caótica. Considere também que a = 0. Se impormos que as condições iniciais dos oscila- 
dores são próximas, porém distintas, teremos que as trajetórias associadas aos vértices irão divergir 
exponencialmente rápido até o ponto em que as diferenças ||xj(í) — Xi(í)|| para todo i, j = 1, ■ ■ ■ ,n 
com i ^ j serão tão grandes quanto o próprio diâmetro do atrator. Em outras palavras, se consi- 
derarmos a rede desacoplada (a = 0) e condições iniciais distintas, então não haverá sincronização 
entre os osciladores. Além disso, se considerarmos agora que a rede está desacoplada mas todos 
os osciladores possuem as mesmas condições inciais, teremos que a variedade de sincronização, sob 
qualquer pequena perturbação, também irá começar a divergir exponencialmente rápido e dessa 
forma levar os vértices a terem comportamentos distintos em suas dinâmicas, ou seja, neste caso, 
a variedade de sincronização perde sua estabilidade sob qualquer perturbação. Não consideramos 
a < O, isso por que neste caso, as soluções de (6.2) podem não ser limitadas, e portanto não faria 
sentido falar de sincronização. Considere por exemplo n = 2, f = OeH = Im- Neste contexto, a 
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diferença xi{t) — X2{t) 



z{t) cumpre 



z{t) = e 



-2at 



z{0) 



de forma que a mesma só irá decrescer exponencialmente em í se a > 0. Portanto, consideramos 
sempre a > 0. 

6.2 Principais Resultados 

A questão da estabilidade da variedade de sincronização N é um dos principais objetivos do pre- 
sente trabalho. Portanto, os nossos principais resultados estão relacionados à mesma. Tais resultados 
são três, e chamá-los-emos de Teorema da Existência Global das Soluções, o qual está relacionado 
ao modelo 6.3, Teorema da Estabilidade da Variedade de Sincronização, que garante a sincronização 
e sua consequente estabilidade e finalmente. Teorema da Persistência da Sincronização, que segue 
praticamente como um corolário do Teorema do Parâmetro Crítico de Acoplamento. Enunciamos 
os mesmos no presente capítulo, porém as suas respectivas provas serão dadas no Capítulo 8. 

A construção da função de Lyapunov, estabelecida no capítulo 3, é usada para garantir a exis- 
tência global das soluções, tanto no contexto das equações desacopladas quanto e principalmente 
para os oscilados em inter ação - estamos nos referindo ao modelo (6.3). Isso porque para tratar de 
sincronização dos osciladores, primeiramente precisamos garantir a existência global das soluções 
e que as mesmas sejam limitadas. Portanto, independentemente da rede considerada, a hipótese 
de que a dinâmica isolada tem uma função de Lyapunov garante que a dinâmica coletiva também 
cumpre a mesma propriedade. Enunciamos então o 

Teorema 16 (da Existência Global das Soluções). Considere o modelo de acoplamento difusivo 



Assuma que a dinâmica isolada tem função de Lyapunov satisfazendo a Suposição 1. Então, inde- 
pendentemente da rede, as soluções do modelo (6.7) certamente entram, em tempo finito, em um 
domínio absorvente Çl, e portanto as mesmas são limitas e existem por todo o tempo futuro. 

A ideia da prova do Teorema da existência global das soluções, que encontra-se na Seção 8.1 do 
Capítulo 8, está baseada na hipótese de que a dinâmica isolada, de todos os osciladores, tem função 
de Lyapunov e consequentemente um domínio absorvente onde as soluções "moram" por todo o 
tempo. Então, tomamos como principal hipótese a Suposição 1. A partir então dessa hipótese 
construiremos uma função de Lyapunov para a dinâmica coletiva. Essa função para a coletividade 
é dada pela soma das funções de Lyapunov para as dinâmicas individuais, mais precisamente, 
consideramos na soma apenas as funções de Lyapunov das dinâmicas cujas trajetórias ainda estão 
fora do domínio absorvente. Note que se todas as trajetórias já estão num domínio absorvente 
então não há nada a se fazer, pois tal domínio é um conjunto compacto, positivamente invariante, e 
portanto a existência global das soluções já estaria automaticamente garantida. Porém, através dessa 
função de Lyapunov para a dinâmica coletiva, conseguimos garantir se se todas as n as trajetórias 
estão fora do domínio absorvente, ou n — 1 ou • • • ou 1, está(ão) fora do domínio absorvente, a tempo 
finito todas as trajetórias, coletivamente, isto é, considerando o acoplamento difusivo, entram num 
domínio absorvente e lá permanecem por todo o tempo futuro. 

Graças a construção da função de Lyapunov para a dinâmica isolada e consequentemente para 
a dinâmica coletiva, sabemos onde moram as soluções. Considerando o modelo de acoplamento 
difusivo (6.3) e sua forma compacta (6.5), condicionamos que a variedade de sincronização é lo- 
calmente atratora se as trajetórias possuem condições iniciais próximas, ou seja, para quaisquer 
i,j = 1, • • • ,n com ||a;j(0) — 3:^(0) || < 5 com í > O pequeno, então Xi{t) — )• Xj{t) \/i,j. Essa con- 
vergência é determinada pela força de acoplamento a entre os osciladores. Além disso, o parâmetro 



n 




(6.7) 
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crítico que produz essa propriedade advém das contribuições da dinâmica isolada dos osciladores, 
de propriedades espectrais da função de acoplamento ií e do segundo autovalor do laplaciano da 
rede. 

Teorema 17 (da Estabilidade da Variedade de Sincronização). Considere o modelo de aco- 
plamento difusivo (6.7). Considere também a existência global das soluções de acordo com o Teorema 
16. Assuma que as condições iniciais estão em uma vizinhança de e considere 

ac = (6.8) 

onde (3 = sup^^Q ||pt£)/(x)P||oo, deforma que a função de acoplamento H do modelo (6.7) é uma 
matriz positiva-definida com a representação espectral H = PDP\ A2 é o segundo autovalor do 
laplaciano da rede e ni é o primeiro autovalor de H. 

Se a > ac, então o estado síncrono é garantido e o mesmo é uniformemente assintoticamente 
estável. 

A prova é apresentada no capítulo 8 Seção 8.2. Abordamos aqui ideia da mesma. Um dos 
principais ingredientes da prova é considerar a análise de uma só equação diferencial em M"™ e 
não de n diferentes equações em M"*. Sob as hipóteses de que as matrizes laplaciana e da função 
de acoplamento são simétricas, é possível realizar uma mudança linear de coordenadas para uma 
nova base ortonormalizada formada pelos autovetores das matrizes anteriormente mencionadas. 
Considerando essa mudança linear de coordenadas, um vetor solução arbitrário é dado por X = 
1 s + U , onde s = s{t) é uma dada solução síncrona fixada e o incremento U é tal X está 
numa vizinhança da variedade de sincronização N. Iremos então conseguir condições para que as 
contribuições dos modos transversais em X convirjam à zero com o tempo, isto é, que os operadores 
de evolução dos modos transversais tenham todos contração uniforme. As nomenclaturas "modos 
normais" e "modos transversais" estão relacionadas com os autovetores da matriz laplaciana da 
rede. O modo normal é exatamente o sub-espaço gerado pelo autovetor 1. Como consideramos que a 
dinâmica individual de cada elemento da rede é essencialmente uma dinâmica não-linear, utilizamos 
a expansão em série de Taylor e analisamos a parte linear do campo de vetores separadamente do 
resto de Taylor. Após conseguir condições para que os operadores de evolução dos modos transversais 
tenham todos contração uniforme, consideramos a influência do resto de Taylor sobre a análise. 
Utilizamos diretamente a Proposição 3 que garante que mesmo considerando o resto de Taylor, 
uma perturbação não-linear, o estado síncrono não sofre alteração. Para conseguir as condições 
mencionadas acerca dos modos transversais, projetamos a equação em blocos nos modos normais e 
transversais afim de analisar apenas as equações diferenciais da mudança linear de coordenadas, que 
por sua vez são equações diferenciais lineares não-autônomas, e utilizamos o Teorema da diagonal 
dominante (11) para finalmente obter uma cota para a força de acoplamento a que por sua vez nos 
dá o critério enunciado. 

O próximo resultado segue praticamente como um corolário do Teorema 17, isso devido ao tipo de 
estabilidade que estamos considerando, isto é, tem-se a propriedade de "robustez" sobre a variedade 
invariante de sincronização global. Assim, devido a forte garantia do Teorema 17 podemos introduzir 
pertubações no modelo considerado e esperar que a estabilidade da variedade de sincronização não 
seja destruída, isso sob certas hipóteses na magnitude da pertubação. 

Vamos trabalhar com um modelo de pertubação que age na função de acoplamento do modelo de 
rede (6.2). Consideramos um modelo de pertubação que age nas funções de acoplamento pois com 
algumas manipulações algébricas é possível utilizar precisamente o Teorema 10, e então estabelecer 
a magnitude que as perturbações podem ter para que o estado síncrono seja persistente. Mais 
precisamente vamos considerar, para cada oscilador 1 < i < n, n funções de perturbação Vij : 

— >• W^, totalizando então n? possíveis funções de pertubação. Consideraremos apenas o caso 
em que cada Vij é um operador linear em IR*", os quais podem ou não ter uma dependência temporal, 
porém a análise é a mesma, e sem perder a generalidade, representamos tais perturbações na forma 
matricial. Considere ainda o modelo de rede (6.2) o qual será escrito na forma do modelo de 
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perturbação considerado como 

n n 

Xi = f{xi) + a ^ AijH{xj - Xí) + ^ AijVij{xj - Xi) 



fixi) + ^Aij{aH + Vij){xj - Xi) 



(6.9) 
(6.10) 



ou ainda, sabendo que Lij = õijgi — Aij tem-se 



f{xi) + '^{ôijgi - Lij){aH + Vij){xj 



n 



= f{Xi)-Y,L^J{aH + V^j){, 

E na forma de blocos, a equação (6.12) é escrita como 



(6.11) 
(6.12) 



X = F{X) - a{L H)X + 



Y,{BiLBj){Dj,-In)®Vij 



(6.13) 



onde Dij £ Mat(R, n) é uma matriz cuja ij-ésima entrada é igual ale todas as demais são iguais O 
(zero). Ai £ Mat(M, n) é a matriz diagonal tal que a i-ésima entrada da digonal é igual ale todas 
as demais são nulas, L G Mat(M, n) é o laplaciano da rede e In G Mat(M, n) é a matriz identidade. 
Maiores detalhes de como a equação (6.13) é obtida estão disponíveis no Apêndice A. 2. 
Temos então o 

Teorema 18 (da Persistência da Sincronização). Considere o modelo de acoplamento difusivo 
(6.2) o qual é reescrito na na forma de perturbação (6.10) e escrito na forma de blocos como descrito 
em (6.13). Considere que a > ac (6.8). Então a variedade de sincronização é persistente se 



sup 

t 



E 



V 



IVíilloo < 211^ 



(6.14) 



onde T] = a\2^i — (3, e m é a dimensão do sistema isolado (6.10). 

A notação i ^ j, significa dizer que o vértice i da rede é vizinho do vértice j, ou então que o 
oscilador i está acoplado com o oscilador j. Antes de enunciar o teorema falamos que é possível 
introduzir ao total funções lineares de perturbação, porém na prática, algumas delas, ou muitas 
delas, não produzirão qualquer influência na estabilidade da variedade de sincronização. Precisa- 
mente, como consideramos que o modelo de perturbação age na função de acoplamento entre os 
osciladores, então só deverá existir perturbação para aqueles osciladores que estão acoplados. 

O principal ingrediente da prova, que estará disponível na Seção 8.3 do Capítulo 8, é a uti- 
lização do Teorema 10. Além disso, a prova de referido teorema segue praticamente os mesmos 
passos da prova do Teorema 17. Ou seja, utilizamos novamente a teoria de equações diferenciais 
lineares não-autônomas abordada no Capítulo 4. O critério sobre a magnitude das pertubações é 
dado então principalmente pelo Teorema 10 e secundariamente pelo Teorema 11. O Teorema da 
persistência da sincronização poderia até mesmo ser enunciado de uma forma mais generalizada, 
isto é, generalizada no que se refere à utilização de uma norma induzida arbitrária substituindo 
(6.14) por uma desigualdade geral, porém se assim o fosse, o seu enunciado envolveria sucessivos e 
dependentes parâmetros os quais aparecem na demonstração do Teorema 17, portanto preferimos 
enunciar o Teorema 18 utilizando apenas a norma infinito. 

Uma das coisas interessantes que a cota (6.14) nos diz é que podemos comparar a "robustez" de 
diferentes redes, isso porque tal cota é dependente e inversamente proporcional a norma infinito da 
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matriz laplaciana da rede. Então, por exemplo, é de se esperar que redes complexas aleatórias sejam 
mais "robustas" do que redes complexas scale-free, pois notoriamente, tem-se que ||-Z^||oo = 2maXí5'j 
onde gt é o grau do vértice i da rede. E como, as redes tipo scale-free tem uma grande quantidade 
de heterogeneidade, ao contrário das redes aleatórias, é de se esperar que essa última seja mais 
resistente à perturbações lineares. Portanto, pelo Teorema 18, podemos comparar a robustez de 
diferente tipos de redes olhando apenas para as suas topologias. De uma forma mais geral, note 
que T] = ri{X2) e que, considerando redes suficientemente grandes, A2 ~ gi (Teorema 14, onde gi é o 
menor grau da rede, tornando o lado direito da desigualdade (6.14) proprorcional à gi/gn e portanto, 
para redes scale-free muito grandes qualquer pequena perturbação poderia destruir a estabilidade 
da variedade de sincronização. 

Queremos enfatizar que a condição (6.14) dada pelo Teorema 18 é apenas uma condição ne- 
cessária, isto é, a mesma não garante que a estabilidade da variedade de sincronização não seja 
resistente a perturbações com magnitudes fora da cota estabelecida. Uma pergunta interessante, 
a qual não iremos tratar nesta dissertação, é: Qual conjunto de perturbações o qual o resultado 
apresentado é justo? Ou seja, quais seriam as perturbações tais que, tendo valores de magnitude 
fora do estabelicido em (6.14), dostroem a estabilidade da variedade de sincronização? 



Capítulo 7 

Ilustrações 



Este capítulo aborda o caso de rede n = 2 que é o exemplo mais trivial de rede. Apesar de ser 
o exemplo mais simplório, será possível extrair a essência dos Teoremas 17 e 18. 

O objetivo então é transformar o problema de estudar a sincronização entre osciladores difu- 
sivamente acoplados num problema de estudar a estabilidade da solução trivial de uma equação 
diferencial linear não-autônoma. Sobre cada oscilador, introduzimos uma cópia da equação (3.1). 
Dessa maneira, todas as variáveis do sistema podem ser escritas da seguinte forma: 

xi = f{xi) +a[H{x2) - H{xi)] , . 

±2 = fix2)+a[H{xi)-H{x2)] ^ ' 

onde a é o parâmetro global de acoplamento. 

7.1 Equação Variacional 

Diferentemente da abordagem feita na prova do Teorema 17, utilizaremos o vetor diferença 
z{t) = xi{t) — X2{t) para estudar a estabilidade da variedade de sincronização, porém, o resultado é 
o mesmo descrito pelo Teorema em questão para n = 2. Utilizando o teorema da existência global 
das soluções (17) garantimos que as soluções xi(t) e X2{t) de (7.1) existem sempre e são limitadas 
valendo-se das hipóteses do referido teorema. Sabendo que as soluções existem globalmente, o 
interesse foca-se na determinação do valor do parâmetro de acoplamento a, tal que a diferença de 
estados, ou seja, xi{t) — X2{t), converge ao vetor nulo. Iremos omitir a dependência temporal das 
variável por uma questão de conforto na notação. Derivando, obtemos 

Z = Xl — X2 

= f{xi) + a\H{x2) - H{xi)\ - f{x2) - a[H{xi) - H{x2)] 
= f{xi)- f{x2)-2aH{z). 

Como por hipótese, devemos ter que ||2;(0)|| ^ 1, podemos realizar a linearização de f{xi) = 
f{x2 + z), baseando-se na expansão em série de Taylor ficando com 

z = f{x2 + z)- f{x2)-2aH{z) 
= Df{x2)z - 2aH{z) + r{z) 

onde Df{x2) é a matriz jacobiana de / em X2- Assim, como por hipótese / é de classe C^, r > 2, 
utlizamos o Princípio da Linearização (Proposição 3), pois é possível garantir, pelo Teorema da 
Taylor com Resto de Lagrange [Apo62] que r{z) = Odlzp), de forma que o resto de Taylor não 
destrói a estabilidade da solução trivial da equação variacional 

z = [Df{x2) - 2aH]z. (7.2) 
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7.2 Parâmetros Críticos de Acoplamento 



Vamos utilizar a Equação (7.2) para estudar a estabilidade da solução trivial z{t) = O, à qual 
representa a variedade de sincronização. Consideraremos a partir deste momento, e por todo o 
restante deste capítulo, que H = 

Queremos estabelecer condições sobre a para que o operador de evolução da Equação (7.2) 
tenha contração uniforme e consequentemente a sua solução trivial seja uniformemente assintotica- 
mente estável. Podemos assim, utilizar um dos critérios para as contrações uniformes abordados. De 
maneira mais específica, vamos utilizar o Teorema da diagonal dominante (11). Portanto, devemos 
ter que 



ou 



a > 



Df{x2h + 



m 

E 



Df{x2\ 



< 



< o 



E p/(^2) 

i=i 



ou ainda, de forma superestimada, e considerando-se que as condições iniciais estão em uma vizi- 
nhança de um domínio absorvente O para a dinâmica /, tomamos 



a > 



1 



sup 

1 < í < m 



e portanto, escrevemos 



onde 



P/2 



/3 = supp/(x) 



(7.3) 



(7.4) 



de forma que o conjunto Í7 é advindo da existência da função de Lyapunov para a dinâmica isolada 
dos osciladores. Como Í7 é um conjunto compacto, pelo teorema de Weierstrass (Teorema 7) segue 
que P sempre existe. Note que devido as superestimações sobre o valor de a, devemos esperar que 
a cota para o Uc que conseguimos seja bem maior do que a realmente necessária para garantir a 
sincronização entre os osciladores. 

Dessa forma, podemos enunciar o resultado, o qual é um caso particular do Teorema 17: 

Teorema 19. Considere o modelo de rede (7.1) com H = Im- Assuma que a dinâmica isolada 
tem uma função de Lyapunov com um domínio absorvente íl. Além disso, assuma que as condições 
iniciais estão em uma vizinhança de Çl e considere 

/3 = sup||D/(x)||oo. 

Seja Uc = PI 2. Então se a > Oc a variedade de sincronização é uniformemente assintoticamente 
estável. 

Foram feitas simulações computacionais utilizando a dinâmica de Lorenz como a dinâmica iso- 
lada dos osciladores. Sabemos que a dinâmica de Lorenz é uma dinâmica caótica, então considerando 
o sistema desacoplado e com condições iniciais distintas, o mesmo não irá sincronizar. Apenas com 
base na intensidade do parâmetro global de acoplamento a, poderemos esperar que a variedade 
de sincronização seja estável. Além disso, como tratado na Subseção 3.1.2, o sistema de Lorenz 
admite uma função de Lyapunov. Dessa forma, tal dinâmica cumpre as hipóteses do Teorema 19, 
considerando as condições iniciais sempre numa vizinhança de Q, neste caso como definido em (3.7). 

Considere o Teorema 19. De posse das principais ferramentas, vamos encontrar o ac para os 
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dois osciladores acoplados. Sendo x = {x,y,z), a matriz jacobiana do sistema de Lorenz é 

Dm = 




e sua norma infinito é dada por ||L'/(x)||oo = max{2o", |r — + 1 + |y| + \x\ + b}, de forma que 
queremos calcular f3 = sup^gQ ||í)/(x)||oo. 

As trajetórias encontram-se dentro do domínio absorvente Q (3.7) que por sua vez é uma região 
limitada por um elipsóide, então o valor máximo atingido por x acontece quando y = O e z = 2r, 
da mesma forma y é máximo quando x = 0ez = 2rezé máximo quando x = y = 0. Assim, 
podemos tomar as seguintes cotas: 

I I ^ 11^ I o I / ''^ 

\x\ < , \y\ < — , e \z — 2r\ < — , 

' " Vb^ ' ' " ^ct(ò- 1) " ^aib-l) 

Mas com relação a z queremos estimar |r — 2;| = |z — r|. Usamos então a desigualdade triangular 
inversa, assim 

rb 

> \z — 2r\ = \ {z — r) — r\ > \ \z — r\ — r\ 



então 



e portanto 



rb 

\\z — r\ — r\ < 



rb rb 

< \z — r\ < — , + r. 



Assim, 



I rb b\fr 6-v/r rb , 1 

/? = max { 2cr, , + r + 1 + , ^ , —^= + ^ + 6 ^ . 

Então, substituindo os valores clássicos dos parâmetros do modelo de Lorenz, ou seja, a = 10, 
r = 28 e ò = 8/3 obtemos 

rb b\fr 

P = , ^ + r + 1 + ^ « 58.22. 

^a{b-\) 

Portanto, o parâmetro crítico de acoplamento é 

= - ~ 29.11. 
2 

A Figura 7.1 mostra a exitência da variedade de sincronização considerando a > Oc- 

Novamente, devido as superestimações para se obter o «c, é de se esperar que o real parâmetro 
crítico para se observar sincronização entre os osciladores seja bem menor do que 29.11 no caso do 
sistema de Lorenz considerando os parâmetros usuais. A Figura 7.3 traz um resultado da simulação 
de ||xi(í) — X2(í)||2 em função de a. Valores de a > 0.5 produzem Yvccít^^o ||xi(í) — X2(í)||2 = 0. 



7.2.1 Utilizando o Critério da Parte Simétrica 

Considere o sistema (7.2). Seja M(x) = Z)/(x) — 2a/3. Fazendo uso do Teorema da Parte 
Simétrica (12), queremos estabelecer condições sobre a tal que a parte simétrica de M(x) seja 
negativa-definida com x = (x, y, z) G íí, para tanto, consideraremos também o Corolário 3. A parte 
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a = 0.1 a = 30 




Xl(í) Xl(í) 

Figura 7.1: Gráficos de Xi(t) x X2{t) com O <t < 100 onde Xi{t) e X2{t) são respectivamente as primeiras 
componentes das trajetórias de xi(í) = {xi{t),yi{t), zi(t)) e X2(í) = {x2{t),y2(t), Z2{t)) no modelo de dois 
osciladores acoplados com dinâmica de Lorenz. Do lado esquerdo com parâmetro de acoplamento igual a 0.1 
e do lado direito com 30. 



simétrica de M(x), aqui denotada por Ps{M{x.)), se lê 



P,(M(x)) = í?M±Mí)!_,„,3, p.5) 



De forma explícita, temos 



-a -2a {r + a-z)/2 y/2 
Ps{M{^)) = \ {r + a - z)/2 -1 - 2a O 

y/2 O -b-2a^ 

Vamos requer que os autovalores de Ps(M{x.)) sejam todos negativos. Note que Ps(M{x)) de- 
pende ainda de y e z, o que significa que o «c é o maior a tal que os autovalores de Ps{M{x)) sejam 
todos negativos com x E Í7. Mas, de posse das cotas \y\ < rb/ \/ aih — 1) e \z — 2r\ < hr / ^ o(}) — 1), 
para cada —rh/^ a(b — 1) < y < rh/ \J a{b — 1) fixado, variamos 2r — br / ^ o(b — 1) < z < 
2r + br/ \/a(b — 1) e a medida que fazemos essa variação, verificamos quais valores de a produzem 
Al < A2 < As < 0. 

Foram feitas simulações computacionais para o procedimento indicado acima resultando, com 
precisão de uma casa decimal, tem-se seguinte parâmetro crítico de acoplamento: 

ac = 13.03. (7.6) 

Agora, utilizando o item 3 do Corolário 3 vamos impor que os determinantes dos menores 
principais de Ps(M(x)) alternam os sinais começando com det(Ps(M(x))i) < O, onde indiciamos 
por P/i(M(x))i, Ps(M(x))2 e Ps(M(x))3 os respectivos menores principais de ordem 1, 2 e 3. A 
ideia é que para cada determinante devemos tomar o maior a do polinómio resultante tal que esse 
polinómio produz valores positivos ou valores negativos dependendo da ordem do menor principal. 
E em seguida, tomamos o maior valor de a dentre os obtidos na etapa anterior, então esse será o 
nosso para este caso. 

Fazendo os cálculos temos 

det(Ps(M(x))i) = -<T - 2a < O 



o que implica a > — (t/2, o que é irrelevante para a análise pois estamos considerando apenas valores 
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positivos. Para o segundo determinante temos 

det(P,(M(x))2) = {-a - 2a)(l - 2a) - ^ > 0. 
Para facilitar a notação faça P2Í<y) = det(Ps(M(x))2). Então 

P2(a) = {-a - 2a)(l - 2a) - ^''^'^"^^ 

= 4a^ + 2{a + l)a + a - l/4(r + s - z)^ 

Note que ^2(0) depende ainda de z. Então o procedimento é o seguinte: verificamos quais as raízes 
de P2{cí) com x G ((3.7)) e tomamos a maior dentre todas as raízes positivas. Isso porque o 
coeficiente de a"^ em ^2(0) é positivo indicando que ^2(0) > O quando a assume valores entre 
{—00, ai) U {aj,+oo), considerando que «j e aj são as raízes de ^2(0) e < Oj. Note ainda que 
P2(o) não possui raízes complexas. Simulações computacionais sugerem que a > 6.5972 produz 
P2{cí) > O com x G íi. De forma mais específica, notamos que a maior raiz de ^2(0) é dada por 

-2{a + 1) + 27(7TiF^4^TFT7^^ 

= ^ , 

onde aj atinge o seu valor máximo quando (r + cr — z)'^ é máximo, implicado que esse valor máximo 
é atingido em z = 2r + rb/y^ a{b — 1). Substituindo os valores de a, ò e r temos 

aj f« 6.5972. 

Por fim, calculando os valores de a que produzem psia) = det(P/j(M(x))3) < O temos 



p^{a) = -8a^+[-46-4(o-+l)]a^+ 



(r + cr — z)^ 
-2a - 2ò(cT + 1) + ^ + ^ 



a-cr6-|--(r+cr-2:) +— ■ 
4 4 



,2 



(7.7) 

E possível verficar que para todo y e z as raízes de ps^a) não assumem valores complexos. Descar- 
tamos as raízes negativas e dentre as positivas verificamos quais produzem ^3(0) < O para cada y 
e z, e dentre essas tomamos a maior. As simulações computacionais produzem o valor a > 7.5546, 
e além disso, que tal valor é atingido quando z = 2r + rb/ sj a{b — 1) e y = —rb/ \J a{b — 1). Assim, 
podemos tomar então 

«c = 7.5546 (7.8) 

que é o maior a que produz todos as desigualdades relacionadas aos menores principais. 

Portanto, temos dois novos parâmetros críticos de acoplamento que são um pouco melhores do 
que aquele apresentado na seção anterior. Porém, a vantagem de utilizar o Teorema da diagonal 
dominante para encontrar o parâmetro crítico de acoplamento é que a mesma análise pode ser 
estendida para redes em geral como pode-se ver na prova do Teorema 17 exibida no Capítulo 8. 

Podemos observar alguns resultados das simulações, os quais são apresentados nas figuras a 
seguir. Para a = 0.1 verificamos que não há sincronização entre os osciladores como mostra a figura 
7.2. Em contraste, para a = 30 verificamos o decaimento exponencial em t. De forma um pouco 
mais precisa temos a relação entre o parâmetro de acoplamento a e a norma euclidiana da diferença 
dos estados xi(í) e X2(í) como mostra a figura 7.3 
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a = 0.1 « = 30 




0.25 



tempo tempo 



Figura 7.2: Do lado esquerdo temos a norma euclidiana da diferença de estados das trajetórias dos oscila- 
dores em função do tempo para a = 0.1 e do lado direito para a = 30. 



a = 0.6 




Figura 7.3: Do lado esquerdo temos agora a = 0.6. Verificamos que o decaimento é mais lento porém o 
sistema tende à sincronização. Do lado direito temos a relação entre a norma euclidiana da diferença dos 
estados contra o parâmetro global de acoplamento. Neste caso, estamos tomando a média temporal da norma 
euclidiana com 1000 < t < 2000 para cada a. Notoriamente, para valores a > 0.5 o sistema irá sincronizar. 



7.3 Perturbações em Dois Osciladores Acoplados 

Da mesma forma que sabemos que o parâmetro crítico de acoplamento pode ser muitas vezes 
maior do que o realmente necessário para se obter sincronização persistente, a magnitude das 
perturbações poderia exceder as cotas estabelecidas pelo Teorema 18 e ainda assim a equação 
perturbada continuar na vizinhança de atração da variedade de sincronização. 

Nesta seção vamos trabalhar com alguns exemplos para o Teorema 18 considerando n = 2, 
H — IfYi e o sistema isolado com a dinâmica de Lorenz. Assim, como ||Zj||oo — ll-^l|2 — 2 e ?7j — 3, 
então o valor estabelecido em (6.14) se resume a 

2 

SUP ^ WVijWoo < I 
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Exemplo 4. Considere o Teorema 18 para n - 
Lorenz e as matrizes de perturbação Vn = V22 



2, H = I3, a dinâmica isolada sendo o sistema de 
- V21 = O G Mat(M, 3) e V12 a matriz constante 



Vi 



12 




Vamos tomar a = 30. Já sabemos que, para o caso que estamos considerando, temos f3 = 58.22, de 
forma que tomamos a = 30 > Oc = (3/2, produzindo assim rj = 30 — 29.11 = 0.89. Então tem-se 



sup 

í 



E 



Wi 



ijWoo 



1^12 1 



que deve satisfazer 2\í^\ < rj/A = 0.2225, implicando que —0.11125 < < 0.11125. Porém, simula- 
ções computacionais mostram que para esta faixa de intervalo, a variedade de sincronização existe 
para a > 0.6 como pode ser observado na Figura 7-4, portanto e certamente para a = 30. Tal figura 
tem como objetivo verificar regiões de sincronização de forma que podemos fazer essa análise através 
de um mapa de cores. Por padrão, a região em preto, no mapa, é a região que verificamos a propri- 
edades citada. As cores no mapa representam a magnitude da média temporal de ||xi(í) — X2(í)||2 
com t entre 1000 e 2000 para cada a e onde xi = {xi,yi,zi) e X2 = {x2,y2,Z2) são as duas 
trajetórias para o sistema de Lorenz, de forma que essa variação acontece continuamente a partir 
de O (preto) até o valor máximo atingido pela norma euclidiana da diferença de estados (amarelo). 



1.8 
1.6 
1.4 
1.2 
B 1 
0.8 
0.6 
0.4 
0.2 
O 




-3 



-2 



o 



Figura 7.4: Mapa de cores para o Exemplo 4. A escala de cores representa ||xi(í) — X2(í)||2. 



Lema 4. Sejam A{t) e B{t) funções matriciais contínuas em um intervalo J C M4. tais que 
\\A(t)\\ < M, ||i?(í)|| < M, e suponha que o operador de evolução da equação x = A{t)x cum- 
pre \\T{t,s)\\ < A:e''(*-"), ?7 < O e > 1. 5e 



t2 



B{t)dt 



< 5 para \ti — Í2I < h, 



então o operador de evolução da equação perturbada y = [A{t) + B{t)]y cumpre 

\\fit,s)\\ < (1 + 5)ke^^'-'^ 
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onde /3 = r/ + 3Mk6 + h'^ In (1 + õ)k. 

Este resultado, cuja prova pode ser encontrada em [Cop78] página 7, mostra que se os coefi- 
cientes das matrizes A(t) e B{t) são limitados, então a estabilidade uniformemente assintótica é 
preservada sobre perturbações "integralmente pequenas". Uma aplicação então que pode ser obser- 
vada é que a estabilidade uniformemente assintótica não é destruída sobre perturbações que oscilam 
suficientemente rápido. 

Exemplo 5. Considere agora as mesmas hipóteses feitas no Exemplo 4 com uma única diferença, 
ou seja, considere agora a matriz de perturbação não- constante 

(Ç cos (wí) ^ cos {ujt) Ç cos /l 1 l\ 

Ç cos {ujt) O O = Ç cos {u}t) 1 O O , 

Çcosjwí) O O / \1 O 0/ 

onde a principio tomamos uo sendo a frequência de oscilação do sistema de Lorenz para os parâmetros 
usuais a = 10, r = 29e6 = 8/3. Verifica-se, através dos dados obtidos pelas simulações que a 
frequência dominante é aproximadamente uj = 4.1888. 

Pelo Teorema 18, para que a magnitude da perturbação não altere a estabilidade da variedade 
de sincronização, faz-se necessário que 

2 

sup \\Vij\\oo = sup IIV12II00 < 7- 

Então, como supj ||Vi2||oo = sup^ 3|Ç|| cos(cjí)| = 3|Ç|, teremos que —7.42 x 10^^ < Ç < 7.42 x 10~^. 
Plotamos, na Figura 7.5, valores um pouco mais realísticos para ^ ea, produzindo novamente regiões 
de sincronização (em preto). 

E perceptível a simetria que encontramos. Tal simetria está relacionada com a dependência tem- 
poral em função do cosseno. Além disso, a inclinação da "língua" que observamos está intimamente 
relacionada com o parâmetro u. 

CO = 4.1888 




-10 -5 O 5 10 



Figura 7.5: Mapa de cores para o Exemplo 5 com uj — 4.1888. 

Simulações computacionais mostram que os mapas de cores para baixas frequências são bastante 
oscilatórios e não simétricos, de forma que não é possível verificar inclinação neste caso. Portanto, 
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se oj é pequeno, isto é, se lo 
constante 



O, então isto equivale a ter uma perturbação V12 com a matriz 



V 




de fato pois lim^j^o cos(a;í) = 1. Iremos considerar então apenas valores "grandes" para lo e apenas 
valores não-negativos visto que a análise seria recíproca pois o cosseno é uma função par. 

A pergunta interessante é: o que acontece então quando u; — >• 00 ? E a resposta é: isso equivale 
a não ter pertuhação, ou seja, tornando V12 = 0. Esta conclusão pode ser verificada pelo Lema 4- 
Note que 



lim 

CJ— í-oo 



Vi2Ít)dt 



lim 3|Ç| / I cos{ujt)\dt 



^.^ smjujt) sgnjujt) ^ ^ 



w 



onde "sgn" é a função sinal. Intuitivamente, quando uj cresce, a perturbação tende a ter uma os- 
cilação muito rápida, fazendo com que a estabilidade da variedade de sincronização, definida por 
xi{t) = X2{t), não seja afetada - de acordo com o Lema 4- 

Além disso, sabemos que para o modelo de dois osciladores acoplados com dinâmica de Lorenz e 
sem perturbação, o real parâmetro de acoplamento que se observa sincronização é aproximadamente 
0.5. Então, a medida que aumenta-se o valor de 00 no modelo de perturbação considerado neste 
exemplo, menor será a inclinação da "língua" no mapa de cores, e portanto quando lo ^ 00 restará 
apenas regiões de não- sincronização abaixo de valores menores que 0.5 para a. Este fato pode ser 
observado na Figura 7.6. 





Figura 7.6: Mapas de cores para diferentes valores de lo, relativo ao Exemplo 5. 



Dessa forma, possuindo dependência em t ou não, nota-se que a variedade de sincronização é 
"bastante" persistente, isso devido a força (magnitude) do parâmetro global de acoplamento a, visto 
que o «c é bem maior do que o real para se observar ao menos a sincronização. Então seguindo 
os mesmos passos é de se esperar que a região de estabilidade uniformemente assintótica pode ser 
maior do que aquela que podemos garantir. 



Capítulo 8 

Prova dos Teoremas Fundamentais 



8.1 Existência Global das Soluções 

Iremos construir uma função de Lyapunov para a dinâmica coletiva de redes de n oscilado- 
res difusivamente acoplados, dada a hipótese de que a dinâmica individual possui uma função de 
Lyapunov. O resultado é o Teorema 16. 

Prova: Omitiremos a dependência temporal de Xj(í) para facilitar a notação. 
Considere que a função V : — )• M dada por 

V{x^) = l/2{xi-a,Q{xi-a)) (8.1) 

onde a £ M™" é fixado e Q é uma matriz positiva-definida, é uma função de Lyapunov para os 
sistemas isolados, i = 1, - ■ ■ ,n com um conjunto absorvente Í7. Considere então o sistema acoplado 
(6.3) que pode ser escrito na forma compacta 

X = F{X)-a{L0H)X = G{X) (8.2) 

onde X = (xi, • • • , X2) e F{X) = (/(xi), • • • , f{xn))- Considere também a função W : M"™ — M 
definida da forma 

H^(X): = |^l/(xí) Vitalquexí gM™\í7 (8.3) 

Afirmação: W é função de Lyapunov para G{X). 

Por construção, W é positiva-definida em relação a f?"™'. Portanto, resta-nos estudar o sinal 
da derivada de W. Temos 3 casos a analisar, porém antes de mais nada, note que considerando a 
dinâmica isolada 

V\xi) = l/2{dxi/dt, Q{xi - a)) + l/2{xi - a, Qdxi/dt) 

= l/2(/(xi), Q{xi - a)) + l/2{xi - o, Qf{xi)) = (/(x,), Q{xi - a)). 

E considerando a dinâmica coletiva tem-se que V'{xi) = {gi{xi),Q{xi — a)) onde gi{xi) = f{xi) — 
aYl]=iLijH{xj). 

Caso 1: Xj G M'" \ Í2 para todo i = 1, • • • , n. Então, temos que 

n 

i=l i j 

= 'Y{f{xi),Q{xi - a)) - aYYLij{H{xj),Qxi) + aYYLij{H{xj),Qa) 

i i j i j 

= Y{f{xi),Q{xi - a)) - g Lijx\QHxj + a^^Y^LijO^QHxj, 

i i j i j 
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porém, note que podemos escrever estas duas últimas somas da forma 

Ujx\QHxj =X\L(E) QH)X e YYl LijO^QHxj = A\L ® QH)X 

i j i j 

onde A = 1 (8) a. Assim, 

A\L QH)X = (10 aY{L ® QH)X 
= {1^ ®a)){L®QH)X 
= il^L) ® {a)QH)X, 

porém, l^L = LI = O, logo Á^[L (g) QH)X = 0. Ficando com 

W'{X) = Y,{f(xi)^ (^i - «)) - «^^(^ ® QH)x. 

i 

Como L é positiva semi-definida e Q e H são positivas-definida então tem-se que L (8) QH é positiva 
semi-definida (Teorema 5). E por hipótese Yliifi^ú^Qi^i ~ ^)) < O visto que Xi £ \ para 
todo i, portanto W'{X) <0. 

Caso 2: Xi G para todo i. Como Q é um conjunto positivamente invariante, então não temos 
nada a fazer pois as trajetórias permanecem em Í7 por todo o tempo futuro. 

Caso 3: Vamos impor, sem perder a generalidade e para facilitar as manipulações algébricas, 
que a = 0. Considere inicialmente que xi £ M™ \ e X2, • • • , x„ G Í7. Assim, 

W'{X) = V'ix^) = {f{xi)-aY,LijH{xj),Qxi) 

j 

= {f{xi),Qxi) - a ^^Lijií(xj),Qxi^ . 
Queremos mostrar nesse ponto que (S^j LijH{xj),Qxi) > 0. Perceba que 



/ Y LijHixj), Qxi )=x\Y LijQHxj = X\L ® QH)X Lijx\QHxj. 

\ j I 3 «=2 3 

Em outras palavras, queremos mostrar que cada parcela da soma X\L (8 QH)X é um número 
não-negativo começando com i = 1, de forma que o argumento utilizado pode ser generalizado. 
Então, sabendo que Lij = ôijki — Aij onde ki é o grau do vértice i, temos 

LijH{xj),Qx}j = {õijki - Aij)H{xj),Qxi^ 

= (^{SijhH{xj),Qxi^ - /^AijH{xj),Qxi'^ 

= ki{H{xi),Qxi) - (^AijH{xj),Qx}j 
= Y {{H{xi),Qxi) - {H{xj),Qxi)) > O 

visto que para j > 1, xj G Í2, então, utilizando a desigualdade de Schwarz [LimlO] é possível mostrar 
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que {H{xi), Qxi) > {H{xj), Qxi) para todo j = 1, • • • , n. Assim, cada parcela da soma 

X^L ® QH)X = Y,Yj Lijx\QHxj 

i j 

é um número positivo, pois basta considerar que o vetor Xj G \ Çl e Xj £ para j = 1, • • • , n, 
j 7^ i. Então o argumento será o mesmo feito para xi. Sendo assim, podemos afirmar que para 
qualquer 1 < k < n, sendo xi, • • • , Xfc S M*" \íl e Xk+i, ■ ■ ■ , x„ G íl teremos que W'{X) < O pois 

k 

W'{X) = Y,V'{x.) 

i=l 

k k 

= ^{f{xi),Qxi) - LijxlQHxj 

i=l 1=1 j 

onde Yli=i Hj Lijx\QHxj > O e por hipótese, Yl'i=i Ui^òi Q^i) < 0. ■ 

8.2 Sincronização Global em Redes 

O resultado a ser mostrado é o Teorema 17. Antes porém considere o seguinte lema cuja prova 
está disponível no Apêndice B.3. 

Lema 5. Considere o sistema (3.1) onde f : M™ — t- M™ é de classe C^, d > 2. Suponha que 
as trajetórias de (3.1) moram num domínio compacto íl C M™. Então, o resto R(s,u) de uma 
linearização de f, em torno de uma curva solução s = s{t) de (3.1), cumpre 

\\R{s,u)\\ < M\\uf 

onde u G é tal que x = s + u está suficientemente próximo de s e M é uma constante uniforme 
em s. 

Ainda antes de iniciar a prova, vamos considerar uma pequena observação sobre o espaço eucli- 
diano M"™". Os espaços M" <^ M™' e M""* são isomorfos e M"™ pode ser decomposto em uma soma 
direta da forma 

^nm = TV e r 

com as projeções 

TTN : M"™ N e vtt : M""" ^ T 

onde TTisi + ttt = Id^nm. Os subespaços N e T C M"'" são determinados pelos mergulhos de M"* e 
j^{7i-i)r)i ijj(^u2Ído pelo laplaciano L da rede. Isto é, tais subespaço são dados por 

iV = {l}(»]R™ e T = {v2,--- ,Vn}(^R"' 
onde {1} representada o subespaço de M" gerado pelo vetor 1 e da mesma forma para {^2, • • • , Vn}- 

Prova (do Teorema 17): Queremos mostrar que a variedade de sincronização N (6.6) é lo- 
calmente atratora. Vamos olhar para as equações do movimento {xi{t) G M"^) não separadamente, 
mas como um único vetor obtido do empilhamento X{t) = (xi(í), • • • ,Xn{t)) G M"™, o que produz 
o modelo escrito na forma de blocos 

X = F{X) - a{L® H)X (8.4) 

com F{X) = (/(xi), • • • , f{xn)), L G Mat(M, n) é o laplaciano da rede. Para facilitar a notação e a 
análise omitiremos a dependência em t de X{t). 
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Como L é simétrica, e por hipótese H também é simétrica, então essas matrizes são ortogonal- 
mente diagonalizáveis (Corolário 1). Considere então a representação L = OMO^ e H = PDP^. 
Da mesma forma, L ^ H tem a decomposição L H = (O (8) P){M D){0 P)^. Então, pelo 
Teorema Espectral, como os autovetores de H formam uma base de M"'" podemos representar 
o vetor X nessa base efetuando a seguinte mudança linear de coordenadas 

X = {0(^ P)Y, 

onde Y = (yi, • • • , yn),yi S W^. Mas, de forma explícita pode-se verificar que 

n 

{O0P)Y = Y,Vi<^Pyi, 

i=l 

onde os Vi^s são os autovetores de L. Estamos considerando que os autovetores de L,vi = 1,V2, - ' ' > 
já estão ortonormalizados. Para facilitar a notação colocamos s = Pyi e escrevemos então o vetor 
X na forma 

n 

X = s + ^^Vi^ Pyi = s + U (8.5) 

i=2 

com U = Y17=2 ® Pyi- 

Introduzimos neste momento a seguinte nomenclatura: chamamos o autovetor 1 de L de modo 
normal e os demais (^2, • • • ,Vn) de modos transversais. Queremos então obter condições para que 
as contribuições dos modos transversais convirjam ao vetor nulo de M""^: 

lim U = 0. (8.6) 

í— >oo 

Isto é, a rede irá sincronizar globalmente pois neste caso tem-se que lim^^oo X = 1 <^ s, o que 
significa 

Xi{t) = X2Ít) = ■■■ = Xn{t) = s{t). 

Reescrevendo a equação (8.4) na nova representação (8.5) e sabendo que [L ® H){1. Ç!) s) = 
{LI (g) Hs) = O temos 

X = F{X) -a{L®H)X 
l®s + Ú = F{l®s + U)-a{LCdH){l®s + U) 
= F{l(^s + U)- a{L H)U. 

Neste ponto, nós assumimos que U é pequeno e realizamos a expansão em série de Taylor de 
F{l<Si s + U) em torno de uma vizinhança da variedade de sincronização A^, onde l<Sis + U pertence 
a essa vizinhança, produzindo 

F{l(g)s + U) = F{1 s) + DF{1 s)U + R{s, U). 

Mas, não é difícil ver que F{1 s) = 1 /(s). Note ainda que DF{X) é uma matriz diagonal em 
blocos onde cada bloco 1 < i < n é a matriz jacobiana Df{xi) o implica que podemos escrever 

DF{l(às) = In(^Df{s). 

Então, ficamos com 

\®s + Ú = l® f{s) + [/„ ^ Df{s)]U - a{L ® H)U + R{s, U). (8.7) 

Considere o operador projeção nos modos normais ttn = 1-1^^0/^, eo seu complemento, o 
operador projeção nos modos transversais ttt = Id^nm — ttat. 

Assim, analisamos a equação (8.7) projetando-a nos modos normais e transversais. Aplicando o 
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operador vr^v tem-se, de um lado 7rAr(l (S" á) + TTj\f(U) = 1 (S" á. E do outro 

TTN (1 «) /(s)) + TTN[In «> Df{s)]U - OTInÍL H)U + TTnR{s, U) = 1 ^ f{s) + TTnRÍs, U), 

e portanto 

l^s = l0f{s)+TTNR{s,U). (8.8) 
Aplicando-se o operador ttt em (8.7) tem-se, de um lado 

itt{1 ® s) + ttt{Ú) = Ú. 

Do outro lado tem-se 

7rT(l«>/(s)) = 0, ■KT[In®Df{s)\U =[In®Df{s)]U e -aTTT[L(àH)U = -a{L®H)U. 
Assim ficamos com 

Ú = [In® Df{s) - a{L (g) H)]U + 7rTfí(s, U). (8.9) 

A equação (8.9) é uma equação diferencial linear não-autônoma com perturbação. Então, como 
queremos conseguir condições para garantir que [/ = O seja uniformemente assintoticamente estável 
em (8.9), faremos o seguinte procedimento: Vamos considerar a respectiva equação homogénea, ou 
seja, 

Ú = [In®Df{s)-a{L®H)]U, (8.10) 

obter uma cota para que o operador de evolução desta equação linear tenha contração uniforme e 
a partir de então, considerar a equação perturbada (8.9) e mostrar que a cota obtida na equação 
linear continua sendo válida mesmo sem desprezar o resto. 

Reescrevendo U na base de T = {v2, ■ ■ ■ , Vn} ^ {wi, • • • , Wn}, isto é, U = Y17=2 ^« PVii onde 
os tí7,'s são os autovetores ortonormalizados de H, ficamos com 



1=2 1=2 

n 

Y,yi®[Df{s)-a\iH\Pyi 



i=2 



OU ainda 



n 



Y^Vid?) {Pm - [Df{s) - a\H]Py,} = 0. 

i=2 

Como todos os autovetores Vi são linearmente independentes, segue que Pyi — [D f (s)—aXiH]Pyi = O 
para todo i, ou simplesmente 

yj = [A{t)-aXjD]yj (8.11) 

com j = 2, • • • , n e A{t) = P'^ Df[s{t))P. Portanto, reduzimos a análise à estudar as equações 
diferenciais dos coeficientes da mudança linear de coordenadas. 

Todas as equações, em (8.11), tem a mesma estrutura diferenciado-se apenas pelo j-ésimo au- 
tovalor de L, \j. Procedemos então neste momento a mesma análise feita na Seção 7.2, utilizando 
o Teorema da diagonal dominante (11). Para que a solução trivial yj = O de cada equação seja 
uniformemente assintoticamente estável devemos impor que 

m 

Aii{t)-aXjDii+ ^ \Aik{t)\<0 Vl<i<m e 2<j<n 

k=l, k=íi 

OU ainda 

^ ^ Aujt) + ELl, \^^k{t)\ ^ ET=1 \^ikit)\ 
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Como os autovalores de L podem ser ordenados da forma O = Ai < A2 < • • • < A„, e, por 
hipótese, H é positiva-definida, então os elementos da diagonal de D podem ser ordenados da forma 
O < ^1 < < ■ ■ ■ < A*m) e as condições iniciais de todas as trajetórias estão numa vizinhança do 
domínio absorvente íi, somos motivamos a tomar o parâmetro crítico de acoplamento 



sup,sn \\P^Df{x)P\\ 
A2AÍ1 



112) 



ou simplesmente Uc = f3/X2^ii, onde /3 = sup^-gQ ||P^D/(x)P||oo- 

Portanto, tomando a > ac, garantimos que os operadores de evolução Tj{t,s) associados aos 
yj's de (8.11) possuem contração uniforme e portanto, o Teorema 11 assegura que 

\\y,m<\\Tj{t,s)y{s)\\ 

<||T,(í,s)||||%-(.)||<fc,e-''(*-)||y,(.)||, 

onde ?7 = aA2/ii — /3. Queremos então, conseguir condições para o operador de evolução T(í, s) de 
U{t) em (8.10) tenha contração uniforme, para tanto, vamos utilizar norma euclidiana, sem perder 
a generalidade, para estimar U{t). Note que ||C/(í)||2 = {U{t),U{t)) e que 

/ n n \ 

{U{t),U{t)) = /^v,^Pyi{t),J2vj^Pyj{t)) 



, j=2 
' n 



Í=2 



^v,^ PTi{t, s)yi{s), Vj ® PTj{t, s)yj{s) 



U=2 



Í=2 



Y.T.''^ ® [yÁs)^T]{t,s)P^]PTiit,s)y^{s) 



v\v, ® yi{s)^T}{t, s)Ti{t, s)yi{s) = ^ 



1=2 



De forma análoga, mostra-se que ||C/(s)||2 = X^"=2 ll^íí^)!!!- 

Como \\U{t)\\l = X;r=2 Ilyi(*)ll2 e \\yi{t)\\l < fcí^e-2''(*"'*)||yí(s)|||, onde os kiS são constantes que 
dependem apenas da dimensão do sistema isolado, então 

n 



1=2 



< 



(.^,) 



-2ri{t-s) 



El 

■í=2 



K2e-2'í(*-^)||[/(s)|| 



onde K = maxifcí. Note que T{t,s)U{s) = U{t). Então \\T{t, s)U{s)\\l < K2e-2'?(*-*)||C/(s)||i, ou 
ainda 



T{t,s) 



U{s) 



Seguindo a definição de norma induzida temos que 



||T(í,s)||2= sup 

II«'IIÍ<1 



r(í,s) 



Uis) 



\\Uis)\\ 



onde w = U{s)/\\U{s)\\l Portanto, ||r(í,s)||2 < fce-'?(*-^). 

Considere a equação (8.9) e sua respectiva equação homogénea (8.10), cujo operador de evolução 
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T{t,s) satisfaz ||T(í,s)|| < Ke~''^^~^\ Como o resto R{s,U) cumpre ||i?(s,í7)|| < M||?7|p em virtude 
do Lema 5, então as hipóteses da Proposição 3 são satisfeitas e portanto a estabilidade não muda, 
isto é, U{t) = O é uniformemente assintoticamente estável. Logo, garantimos que a variedade de 
sincronização é uniformemente assintoticamente estável. ■ 



8.3 Perturbações na Função de Acoplamento 



Provamos o Teorema 18. 

Porva: Considere o modelo de perturbação (6.10). Como os operadores de perturbação Vij : M™ — í- 
M™, são operadores lineares, consideramos, sem perder a generalidade, que para todo i e j, Vij é 
uma matriz de dimensão m. Então, na forma de blocos o modelo (6.10) é escrito como 



X = F{X) - a{L(S) H)X + 



^{BiLBj){Dji-In)®Vi 



X, 



5.13) 



onde X = (xi,--- ,Xn), F{X) = {f{xi),--- ,f(xn)), L é o laplaciano da rede, Bi £ Mat(M, n) é 
uma matriz diagonal cuja z-ésima entrada da diagonal é igual ale todas as demais entradas são 
O (zero) e Dij G Mat(M, n) é uma matriz cuja ij-ésima entrada é igual ale todas as demais são 
nulas. 

Realizamos uma análise similar à feita na prova do Teorema 17, isto é, utilizamos a representação 
X = l0s-\-U em (8.13) e analisamos a parte linear do campo de vetores F com a influência do resto 
de Taylor, de forma que procedemos a linearização de F(l (8) s + C/) em torno de 1 (8> s considerando 
que 1 fSi s + U pertence a uma vizinhança de 1 ig) s. Assim 

liE)s + Ú = F{l<E)s) + DF{1 ® s)U -a{L<S)H){l(g)s + U)+ 



+ 



Y,{BiLB,){D,i-In)®Vi, 



s + U) + R{s,U) 



ou 



l®s + Ú = 1® f{s) + [In ® Df{s)\U- a{L H)U + 



3^ 



Vi 



U + R{s,U). 
(8.14) 

visto que F(l s) = f{s), DF{1 ®s) = In® Df{s), a{L ®H){l®s) = a{Ll H{s)) = O e 



{B,LBj){D 



Y,{B,LBj){D 



V,, 



{l®s) = 



pois {BiLBj){Dji — In)! = {BiLBj)Djil — {BiLBj)l onde {BiLBj)Dji é uma matriz diagonal onde 
a i-ésima entrada da diagonal é o número Lij e todas as demais entradas são nulas e portanto o 
vetor {BiLBj)Dji\ tem na sua i-ésima entrada o número Lij e todas as demais entradas nulas. Por 
outro lado, o vetor {BiLBj)l é exatamente o vetor que tem o número Lij na sua i-ésima entrada e 
todas as demais nulas. Portanto (BiLBj)Djil — {BiLBj)l = 0. 

Novamente, fazemos as projeções de (8.14) nos modos normais e tranversais aplicando os opera- 
dores tttv e ttt como descritos na Seção 8.2. Aplicando ttm ficamos com l(8)é = 1® /(s) -|-7rAri?(s, U). 
Aplicando ttt ficamos com 



Ú=[ln0 Df{s) - a{L (g) H)]U + 



Y,{BiLBj){Dji-In)®V 



U + ttnR{s,U). 



5.15) 
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Porém, perceba que se olharmos apenas para a equação não-per turbada 

Ú=[InCS> Df{s) - a{L » H)]U 



(8.16) 



já sabemos, pelo Teorema 17, que U{t) cumpre limt^oo U (t) = O quando a > /3/(A2/^i). Então, 
consideramos a equação (8.16) com a perturbação linear Y17j=ii^i^^j)(^ji ~ ^n) ® Vij U, uti- 
lizamos diretamente o Teorema 10 para determinar a magnitude que a referida perturbação pode 
ter. Assim ■^^{BiLBj){Dji — In) ® Vij deve cumprir. 



sup 

t 



Y,{BiLBj){D,i-In)CSV,, 



< 



V 



onde rj = A2/X1 — /3 e n = max2<í<n ki (ver página 52 para mais detalhes). Mas 

n 

< sup \mLBj){Dji - In) » Vi 



sup 

t 



Y,{BiLBj){Dji-In)®Vi, 



sup^ \mLB,){Dj,-in)\mj\ 



< sup ||L||(p,-,|| + 11^11) ^ \\Vij\\, 
onde ||i?jLi?j|| < ||L|| ||, \\Bi\\ = 1, ||-Djj|| = 1, e = 1, isto é, podemos impor que 



sup ^ II Vij II < 



V 



2k L 



ou melhor, sup ll^íill < 



T] 



2k\\L\ 



il7) 



pois, olhando para o modelo de perturbação (8.13) verifica-se que quando i não é vizinho de j então 
BiLBj = O € Mat(]R, n) de forma que a equação (8.13) pode ser melhor escrita como 



X = F{X) - aiL<S)H)X + 



^ iBiLBj){D,i-In)ç^Vi 



X 



e ainda, as contas realizadas não sofrem alterações. Pelo Teorema 11, tem-se que k = 1 se usarmos 
a norma || • ||oo. Considerando agora a influência do resto na equação (8.15), utilizamos novamente o 
Lema 5 e a Proposição 3, de forma que garantimos que a estabilidade da variedade de sincronização 
não é destruída. ■ 



Capítulo 9 

Conclusões 



9.1 Considerações Finais 

Ao se abordar o problema da sincronização em redes complexas utilizamos a teoria de contra- 
ções uniformes em equações diferenciais lineares não-autônomas, caracterizando assim uma nova 
abordagem sobre um tema tão importante e que está enraizado nos principais ramos das ciências 
naturais, exatas, sociais, e enfim na vida humana. 

Observa-se que o estado síncrono global dos osciladores difusivamente acoplados é uma variedade 
invariante pelas equações do movimento ao considerar o modelo utilizado (6.3). Tal estado é atingido 
e garantido ser uniformemente assintoticamente estável devido à força do parâmetro global de 
acoplamento q, de forma que o seu valor crítico é determinado unicamente pelas contribuições 
da dinâmica intrínseca dos elementos, de propriedades espectrais da matriz de acoplamento e do 
segundo autovalor do laplaciano da rede. 

O parâmetro crítico de acoplamento, com relação à sincronização, é observacionalmente, e como 
já era de se esperar, devido as várias superestimações sobre o mesmo, muito maior do que o realmente 
necessário para que tal estado ocorra. Porém, não se garante que durante o intervalo entre o real 
Qíc e /3/(A2/Ui) o estado síncrono global seja uniformemente assintoticamente estável - a variedade 
invariante (6.6) é localmente atratora se a > /3/(A2^i) (Teorema 17). E também possível verificar 
esse resultado através das simulações realizadas em dois osciladores acoplados. 

Da mesma forma que verifica-se a superestimação de Oc, a vizinhança de atração de N pode ser 
bem maior do que aquela à qual podemos garantir. Além disso, a cota (6.14) para a magnitude da 
perturbação pode apresentar algumas falhas, uma delas está relacionada à utilização da desigual- 
dade triangular para fazer superestimações. Além disso, a condição (6.14) é apenas uma condição 
necessária, isto é, a mesma não garante que a estabilidade da variedade de sincronização não seja 
resistente a perturbações com magnitudes fora da cota estabelecida. 

9.2 Sugestões para Pesquisas Futuras 

A principal sugestão seria a abordagem da persistência da sincronização em redes de sistemas 
não-idênticos. Este caso já foi estudado por Pereira e colaboradores em [PEB''"13]. A variedade 
invariante N não existe neste caso, mas pode-se utilizar resultados da teoria de equações diferen- 
ciais não-autônomas para estabelecer a persistência do comportamento coletivo sob perturbações 
intrínsecas na dinâmica individual dos elementos. 

Esta situação é mais desafiadora por várias razões. Por exemplo, como mencionado, o estado 
completamente sincronizado não é um estado invariante do sistema, e portanto os elementos po- 
dem não sincronizar mesmo para acoplamentos fortes. Pode-se considerar que a não-identidade na 
dinâmica dos osciladores é na verdade uma quase identidade, ou seja, seria o mesmo que considerar 
que o campo do sistema isolado está sujeito a perturbações e então mapeamos 

fiXi) ^ f{Xi) + Wi{Xi), 
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onde maxj ||ií;j(xi)|| < 7 seria a máxima magnitude das perturbações. Mesmo se todos osciladores 

tem um estado inicial igual a perturbação Wi, eles podem ser levados a estados completamente 
distintos. Uma questão em aberto na literatura é se as perturbações do campo de vetores implicam 
na não existência de soluções síncronas. 



Apêndice A 

Álgebra Linear 



A.l Menores Principais 

Nós podemos aferir a definitude^ de uma matriz simétrica calculando os autovalores da mesma. 
Porém, outra forma consiste em usar os menores principais. 

Definição 31. Seja A = [Aij] G Mat(]R,n). Para cada inteiro 1 < k < n, o k-ésimo menor 
principal de A, denotado por D^, é o determinante da matriz 



A, 



Au Ai2 
A21 A22 



A2k 



_Aki Ak2 ■ ■ ■ Akk_ 
a qual se obtém deletando-se as linhas e colunas de ordem k + 1 até n. 
Teorema 20. Seja A G Mat(M, n) uma matriz simétrica e 1 < k < n. Então 

1. A é positiva- definida se e somente se Dk > O para todo k. 

2. A é positiva semi-definida se e somente se Dk > O para todo k. 

3. A é negativa- definida se e somente se { — \)^Dk < O para todo k. 

4- A é negativa semi-definida se e somente se { — \)^Dk < O para todo k. 
Prova: Ver [Zha99], p. 200. 

Corolário 3. Seja A G Mat(M, n) uma matriz simétrica. As seguintes afirmações são equivalentes: 

1. A é negativa-definida. 

2. Todos os autovalores de A são negativos. 

3. Os menores principais de A alternam o sinal começando com Di = det{Ai) < 0. 



^Consiste da classe dos conceitos das matrizes positiva-definida, positiva semi-definida, negativa-definida e negativa 
semi-definida. 
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A. 2 Obtenção do Modelo de Perturbação em Blocos 

Esta secção trata sobre a obtenção do modelo de perturbação na função de acoplamento o qual, 
quando escrito na forma de blocos é resumido na equação (6.13). 
Considere a equação (6.12) que pode ser escrita como 



fixi 



a 



a 



^ LijH{xj) - a ^ LijH{xi) - ^ LijVijXj + ^ LijVij 
j j j j 

LijH{xj) — ^ LijVijXj + ^ LijVij Xi. 



(A.l) 
(A.2) 



O objetivo então é empilhar os n vetores de M"^ de cada parcela de (A.2). Quando do empilhamento 
da percela LijH{xj), isto é, o vetor (Ylj LijH{xj)^ ■ ■ ■ , Y2j LnjH{xj)^ £ M"™", já sabemos que 
o mesmo pode ser escrito como (L iS" H)X onde L é o laplaciano da rede e X = (xi, ■ ■ ■ ,Xn), 
Xi £ M"*, i = 1, - ■ ■ , n. Do empilhamento dos vetores Y2j ^ij^ij^j, pode-se observar que o mesmo é 
obtido através da multiplicação 



/ -í^ll^ll -^12^12 
L21V21 -Ív22^22 

\LnlVn2 Ln2Vn2 



LlnVin\ 
-í^2nV2n 



X2 



Mas a matriz acima descrita pode ser escrita como j=i Dij0LijVij ou ainda como BiLBj® 
Vij onde Dij G Mat(M, n) é a matriz cuja ij-ésima entrada é igual ale todas as demais são nulas 
e Bi € Mat(M, n) é uma matriz diagonal tal que a i-ésima entrada da diagonal é igual ale todas 
as demais são nulas. Note que o produto BiLBj resulta numa matriz onde o ij-ésimo elemento é 
igual a Lij e todos os demais elementos são nulos. 

De forma semelhante, o empilhamento das parcelas LijVijXi pode ser escrito pela multipli- 
cação 



o 



o 



o 
o 



X2 



V 



onde O G Mat(M, m) é a matriz nula. Esta última matriz pode ser escrita como ^^Li -^i' 



Yij=iDiiLij (g) Vij, ou ainda como j=iiBiLBj)Dji (g) Vj. Note que o papel 



de Dji é de trazer o único elemento não nulo, provavelmente, de BiLBj para a diagonal de sua 
respectiva linha. Assim, de posse das parcelas do bloco temos o modelo 



X = F{X) - a{L(S) H)X - 
= F{X)-a{L0H)X + 



BiLBj (g) Vi 



X + 



{BiLBj)Dji ® Vij 



X 



Y{B,LBj){Dj,-In)(^Vi 



X. 



Apêndice B 

Equações Diferenciais Ordinárias 



B.l Existência, Unicidade e Extensão 



Considere o sistema de equações diferenciais autónomas 



X = f{x) 



(B.l) 



onde f : U C — )• M™ é de classe C , ou seja, suas derivadas parciais primeiras existem e 
são contínuas. Abordamos apenas o caso autónomo, porém os resultados aqui apresentados são 
extensíveis ao caso não-autónomo, isto é / = f{t,x). 

Uma solução de (B.l) é o caminho x : J — )• definido para algum intervalo J C M tal que 
x{t) = f{x{t)) \f t £ J. Um valor inicial para a solução x : J — )• M"* é uma especificação da forma 
x(to) = xq onde íq £ J e xq £ W^. Podemos considerar íq = 0. 

Uma equação diferencial não- linear pode ter várias soluções que satisfazem o mesmo valor inicial. 
Além da questão da unicidade, também temos a questão da existência. Então, afim de estabelecer 
condições para que o sistema (B.l), sob uma dada uma condição inicial, possua uma única solução, 
consideramos o 

Teorema 21 (da Existência e Unicidade). Considere o problema de valor inicial 



onde Xq S M™. Suponha que f : U C M™ — )• é de classe C . Então, exite uma única solução 
para o referido problema de valor inicial. De forma mais precisa, existe c > O e uma única solução 



satisfazendo a condição x{tí) = xq. 

Prova: Ver [HSD04], p. 385. 

O resultado apresentado acima garante a existência e unicidade da solução, porém apenas lo- 
calmente. O resultado a seguir dá as condições necessárias para que a solução possa ser estendida, 
ou seja, para que a solução esteja definida não apenas para uma vizinhança do tempo inicial, mas 
por todo o tempo futuro. 

Teorema 22 (Extensão). Seja C um subconjunto compacto de U C e seja f : U C M™ — )• W^, 
f £ C^. Seja xq £ C e suponha que toda solução da forma x : [0,b] ^ U com x{0) = xq encontra-se 
inteiramente em C. Então existe uma solução x : [0,oo) — í- U satisfazendo x{0) = xq e x{t) £ C 
para todo t > O, e portanto a solução está definida por todo o tempo futuro. 



X = f{x) 



x(0) = Xq 



X : (-c,c 



;) 



Prova: Ver [HSD04], p. 397. 
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B.2 Desigualdade de Grõnwall e Variação dos Parâmetros 

O lema que segue é de fundamental importância para obter várias estimativas. A sua prova pode 
ser encontrada em [AM05], página 29. 

Lema 6 (de Grõnwall). Sejam uj{t), k{t) : J C M — t- M+ funções contínuas com oj{t), k{t) > O para 
todo t & J e C uma constante não-negativa. Se 



ijj{t) < C + / k{u)uj{u)du Vs,i 

uj{t)<Cexp( í k{u)dv\ (B.2) 



J s 

então 
para todo t £ J. 

A seguinte proposição também é utilizada em algumas demonstrações, justificando portanto a 
sua abordagem. 

Proposição 8 (Variação dos Parâmetros). Sejam A : M — )• Mat(M, m) e 51 : M — )• M™" funções 
contínuas. Considere a equação não-homogênea 

X = A{t)x + g{t). 

A solução de tal equação com a condição inicial a{s) = xq é dada por 

x{t) = T{t, s)x{s) + / T{t,u)g{u)du, 

J s 

onde T{t, s) é o operador de evolução da equação homogénea correspondente. 
Prova: Ver [Tcsl2], p. 81. 



B.3 Prova do Lema 5 



Prova: Tomando x = s + u suficientemente próximo de s, seguimos a expansão em série de Taylor 
de f{s + u) até ordem 1. Vamos utilizar a representação / = (/i, • • • , fm) onde fk '■ M*" — )• M, A; = 
1, • • • , m. Assim, na expansão, cada componente de / é da forma /^(s + n) = fk{s) + Dfk{s)u + 
rfc(s, u), isto é, 

f{s + u) = (/i(s),-- - ,frn{s)) + {Dfi{s)u,--- , D fm{s)u) + {ri{s, u) , ■ ■ ■ ,r„(s,M)). 

Pelo Teorema de Taylor com resto de Lagrange [Apo62] tem-se que ^^(s, u) = {l/2)u^ D"^ fk{s + 6u)u, 
para algum O < 6^ < 1, Vfc. Assim, ||rfc(s,tí)|| < {l/2)\\D'^fkis + 9u)\\\\u\\'^. Como fk G C"^,d > 2, 
então D"^ fk{s + 6u) é contínua e s + 9u Çi Vi (compacto). Então, pelo Teorema de Weierstrass 
existe Mg tal que ||D^/fc(s + 6u)\\ < Ms para cada s fixado. Assim, Mg : M™" — )• M+ é uma função 
contínua que toma valores num compacto, logo, existe M, uniforme em s, tal que ||rfc(s,'u)|| < 
M||ti|p, VA:. Podemos utilizar db norma || * ||oo; SGin perder a generajlid.ajd.e, para estimar Ri^s^ii) — 
(ri(s,u),-- - ,rm{s,u)). Então 

||-R(s,íi)||oo = max ||rfc(s,ti)|| < M||u||^. 
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